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INTRODUCTION 
1 . Gener a i 

The promise of filaunentary composite materials, whose de- 
velopment may be considered as entering its second generation, 
cont.’nues to generate intense interest and applications activ- 
ity. Such interest and activ ty arc well founded, since they 
are based on the possibility of using relatively brittle mate- 
rials with high modulus, high strength, but low density in 
composites with good durability and high tolerance to damage 
and which, when they do fai] , do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantia ' , y improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the mid 1960 's. Rather limited applications to 
prxiAary aircraft structure have been made, however, mainly in 
a material-substitution mode on military aircraft, except for 
a few experiments now being conducted on large passenger air- 
planes . 

To fulfill the promise of composite materials complete- 
ly requires a strong technology base. NASA and AFOSR recog- 
nize the present state of the art to be such that to fully 
exploit composites in sophisticated aerospace structures, the 
technology base must be improved. This, in turn, calls for 
expanding fundamental knowledge and the means by which it can 
be successfully applied in design and manufacture. 
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As technology of composite materials and structures moves 
toward fuller ^option into aerospace structures, some of the 
problems of an earlier era are being solved, others which 
seemed important are being put into perspective as relatively 
minor, and still others unanticipated or put aside are emerg- 
ing as of high priority. The purpose of the RPl program as 
funded by NASA and AFOSR has been to develop critical advanced 
technology in the areas of physical properties, structural 
concepts and analysis, manufacturing, reliability and life 
prediction. 


2. The Rt^I Program 

Our aooroach to accomol ishira these aoals is throuqh an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing to the logical 
conclusion of such activities; namely research into the com- 
posite structure's long term integrity under conditions perti- 
nent to service use. Inherent in the RPI program is the moti- 
vation which basic research into the structural aspects pro- 
vides for research at the materials level, and vice versa. 

Second, interactions among faculty contributing to program 
objectives - which is a group wider than that supported under 
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the project - ia on a day to day basis* regardless of organi- 
zational lines. Program management is largely at the working 
level, and administrative, scientific and technical decisions 
are made, for the most part, independent of considerations 
normally associated with academic departments. Involvement 
of this kind includes - depending on the flow of the research - 
faculty, staff and students from chemistry, chemical engineer- 
ing, civil engineering, materials engineering, and the depart- 
!Monf of mechanical engineering, aeronautical engineering and 
mechanics . 

Botfi of these characteristics of the NASA/AFOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

The [)rogram objective of advancing fundamental knowledge 
in the areas of physical properties, structural concepts and 
analysis, manufacturing, reliability and life prediction, con- 
tinues as our program evolves, but changes in direction will 
be institute! in the next six months recognizing the evolution 
of the use of composites in the industry and among certain 
government agencies such as FAA, DOD and NASA itself. 

Research on new structural design concepts at RPI will 
move away from redesign of actual aircraft components, such 
as the 727 elevator actuator attachment rib and the L-1011 
engine drag strut (formerly known under the acronym CAPCOMP 
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for Composite Aircraft Program Comp onent) , toward those spe- 
cific subjects which these efforts have revealed as areas of 
critical technological need. 

Currently, such specific subjects include the following: 

(a) attempts to quantify St. Venants' principle 
for typical composite materials and member cross 
sections. This is needed to determine hv'^w far 
from points t ;oncentrateu load introduction 
(such as occur in mechanical joints) structural 
cross sections must vary before constant cress 
sections can be efficiently used. 

(b) research into scaling laws for composite 
structures. Here at least two aspects introduce 
uncertainty; first, the effect of constituent mate- 
rial property differences with size as related to 
cure processing exothermic reactions, and second, 
the effect of the ratio of fiber diameter to the 
smallest dimension in the scaled part. The impor- 
tance of advancing understanding in this ar a is in 
the usefulness and efficiency of experiments in the 

cvelopment of large, complex structures such as 
arc the stock-in-trade of the aerospace industry. 
Developments using scaled test specimens would re- 
duce material costs; ease handling, manufacturing 
equipment and tost equipment requirements; and 
generally shorten development schedules, 

(c) lug design using composites. In this area, 
the favorable strength to weight ratio of compos- 
ite materials tends to be offset by their less 
than favorable strength to volume ratio. The com- 
plex load-pathi, furthermore, tend to frustrate 
attempts to use the unidirectional advantages of 
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"tailored materials". Considerable promise remains, 
however, in the possibility of specially "tailoring" 
material properties both azimuthally and radially 
from the points of load introduction. Much remains 
to be done of a fundamental nature in this area both 
theoretically and experimentally. 

Overall program emphasis, in future months - on all as- 
pects - will move from specific applications to more basic, 
long-term research in the following categories: (a) constitu- 

ent materials, (b) composite materials, (c) generic structural 
elements, (d) processing science technology and (e) maintain- 
ing long-term structural integrity. Those computer software 
developments will b»' pursued which both support Rensselaer 
projects in composite materials and structures research n the 
areas listed above and which also represent research with the 
potential of widely useful results in their own right. 

The sailplane project (called CAPGLIDE for Composite Air- 
craft Program Glide r) involving design, fabrication and testing 
of an advanced, full-scale, flight-worthy, all-composite air- 
craft will be continued as a highly productive way to conduct 
research on processing science and technology for low-cost 
aerospace structures. The list of innovative designs and pro- 
cesses resulting from that research to date includes: 

a) A simplified fabrication process without molds 
for accurate, smooth surface airfoils in wings 
and tail surfaces up to roughly twenty feet 
semi-span, 

b) Single sided, composite-foam "sandwich" struc- 
tures which can be successfully joined to 
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sub-structures to form ultra-light weight sur- 
faces capable of carrying useful distributed 
loads , 

c) Composite hinge configurations that are integ- 
rally bonded into airfoil contours to minimize 
pee 1 -strength problems, 

d) Proper mixtures of hardeners and experimental, 
especially low viscosity resins to achieve 
both minimum-depth bond and surfaces for 
ultra-light weight and acceptable strength, 
ami 

e) A i-ochnigue for fabricating smooth-surfaced 
graphite/epoxy tubes of relatively small dia- 
meter, witliout the use of release agents or 
external molds, through the use of teflon 
rods as mandrels. 

The impetus toward innovation in processing science will 
be preserved as the project evolves, by undertaking continu- 
ally more challenging structural functions in sailplane de- 
sign. Further, research on problems within the disciplined 
context of actual manned flight testing can be expected to 
contribute to fundamental knowledge as regards monitoring and 
maintaining structural integrity of functioning, low-cost aero- 
space structures. 

The changes in program direction discussed above will not 
begin to take effect until tlic six-month period immediately 
ahead. This report, however, presents progress achieved in 
all elements of the program under a format better suited to 
the new program orientation, as reflective of converting to 
operation in the new mode. 
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Accordinqly, proqram element hoadinqs In the followinq 

paqes will be 

Constituent Materials, 

Composite Materials, 

Generic Structural Elements, 

Processinq Science and Technology and 
Maintaining Structural Integrity. 

J . Technical Interchange 

Tc^chnical meetings, c^n- and off-campus, provide important 
opportunities for intc^rchange of technical information. Be- 
c'ause of the lanje number of c'ornposites meetings, a c'entral 
catalog with all u, -coming meetings is being maintained and 
distributc'd j'-'r i od i Cti 1 1 y . In this way we help assure that a 
Ronssc'laer facul ty/st .»! f member can participate in important 
mec’tings. The’ caletular for this reporting period is shewn in 
•j’able 1-1 , Mec’tings attc»ndc'd by RPT ccamposites proqram facul- 
ty/staff/students during the' re'porting period are shown in 
Table 1-?. Some meetings [particularly relevant to composites, 
lu'ld en-e'ampus with special s[pe'ake'r s , are listed in Table 1-3. 
A list of ccPUH’cPS i t e'-rc>late'd visits to rele'vant organizations 
by RPI f ae'u I ty/sta f f /st ude'nt s , with the purpose of each visit 
outlined, vs {irese'nte'd in Table I -4. 

Two coivtinuing e'ducat ion sjpecial courses, outc^rowths of 
the composite' mate'rials and structurevs program, were? pre- 
sentc'd for graduate' enqinee'rs in industry and government dur- 
ing the' month evf July. They wore "Composite Materials and 
Structure's" and "The Finite Element Approach in Structural 
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CALENDAR OF COMPOSITES-RELATED MEETINGS 
(April 30, 1981 throuijh Scpteir.ber 30, 1981 ) 


5/12-13 

5/12-14 

5/17-20 

5/18-21 

5/18-21 

8/2-4 

6/5-7 

6/15-17 

6/15-19 

6/21-26 

6/22 

6/22-24 

6/22-26 

8/17-21 

8/26-27 

9/8-9 

9/16-18 

9/20-23 

9/22-25 


6th ASTM Conference on Composite Materials: Testing 

and Design, Phoenix, AZ . "Sponsored by ASTM Committee 
D-3 0." 

50th Anniversary Meeting and Technical Display, Long 
tU'ach, California. "Sponsored by AIAA." 

37th Annual Forum and Technical Display, New Orleans, 
LA. "Sponsored by AHS . " 

international Materials Policy and Technical Consider- 
ations for Future Automotive Design, Dearborn, MI. 

International Symposium on Mechanical Behavior of 
Structural Media, Ottawa, Canada. 

Colloquium on Advanced Mechanics of Reinforced Con- 
crete, Delft, Netherlands. "Sponsored by ASCE." 

Spring Meeting (Electro Dynamics in Solid Mechanics, 
Nonlinear Dynamics) , Alexandria, VA. "Sponfored by 
S I AM . " 

3rd Symposium on Dynamics and Control of Large Flex- 
ible Spacecraft, Blacl^sburq, VA. "Sponsored by VPI 
and SU/AIAA." 

U.S.A.-ltaly Joint Symposium on Composite Materials - 
The Role of Polymeric Matrix in Their Processing and 
Structural Properties, Capri, Italy. 

Pressure Vessels and Piping Conference, Denver, CO. 
"Sponsored by ASME." 

ASEE National Meeting, Los Angeles, CA. 

Joint ASMK/ASCE Mechanics Conference, Boulder, CO. 

15th Biennial Conference on Carbon, Philadelphia, PA. 

6th International Conference on Structural Mechanics 
in Reactor Technology, Paris, France. 

National American Clu'mical Society Meeting, N.Y.C., 

NY . 

7th Europeari Rotorcj-aft and Powered Lift Aircraft 
Forum, Garmisch, Ger.nany. 

International Conference on Composite Structures, 
Paisley, Scotland. 

4th Failure Prevention and Reliability Conference 
(Joint ASME Vibrations, Design Automation, F-P , R 
Conference) , Hartford , CT. 

Euromech 147, Symposium on Damage Mechanics, Cachan, 
Franco. 
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i"?. 

C OMINGS I TES- RE LATl'U) T ECHNICAL MEETINGS ATTENDED OFF-CAM PUS 
(April JO, 1981 through September 30, 1981) 


5/12-13 Oth ASTM Conference on Composite Materials: Testing 

and Design (Prof. Krcmpl) , Phoenix, AZ 

ASTM E9.01 Committee Meeting: Professor Krempl 

presented the paper, "Time Dependent Effects 
During Cyclic Loading of Structural Steels at 
Room Temperature" . 

ASTM E9.08 Committee Meeting: Professor Krempl 

presented the paper, "Preliminary Experience 
witli the Instron Biaxial Extonsometer " . 

5/13-16 50th Anniversary ATAA Meeting and Technical Display 
(Prof. Loewy) , Long Beach, CA 

6/15-19 U . S . A . - 1 1 a 1 y Joint Symposium on Composite Materials - 
The Role of the Polymeric Matrix in Their Processing 
and Structural Properties (S. Sternstein) , Capri, 

I t a I y 

Piofessor Sternstein, Invited Lecturer, pre- 
sented the paper, "Properties of iji situ Epoxy 
Ma t r i ('es" • 


ti/22 ASEE National Meeting: Education & Industry - .A Joint 

Ende.woi' {Vvo( . Shephas'd) , Los Angeles, CA 

Piotessoi Shephard presented the paper, "The 
University's Role in Meeting Industry's Role 
in CADCAM Needs". 

6/22-24 Joint ASMl'./ASCE Mechanics Conference (Profs. Hoff and 
Krempl and student C. Muser), Boulder, CO 

Profe.ssor Hoff presented the paper, "Stress 
Concentrat ions in Cyl indr ical ly Orthotropic 
Composite Plates with a Circular Hole". 

Professor Krempl presented the paper, "The 
Rate (Time) -Dependent Behavior of Ti-7Al’- 
2Cb-lTa Titanium Alloy at Room Temperature 
under Quasi-static Monotonic and Cyclic 
Load i ng " . 

6/22-26 15th Biennial Conference on Carbon (Prof. Diefendorf), 
Ph i ladelphic , PA 

Professor Diefendorf presented the paper, 

"Tlie Chemical Vapor Deposition of Carbon 
Ma trices" . 
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TABLE 1-2 (continued) 


8/17-21 fcth International Conference on Structural Mechanics 
in Reactor Technology (Prof. Krempl) , Paris, France 

8/26-27 National American Chemical Society Meeting (Profs. 
Diefendorf and Sternstein) , N.Y.C., NY 

Professor Diefendorf presented the paper, 
"Consequences of Residual Stress on Carbon 
Fiber Performance". 

Professor Sternstein presented the paper, 
"Viscoelastic Characterization of Neat 
Rosins and Composites" . 

9/20-23 4th Failure Prevention and Reliability Conference 

(Joint ASME Vibrations, Design Automation, F-P, R Con- 
ference) (Prof. Brunelle and L. VJ. Chen and L. Y. 
Chen), Hartford, CT 

Professor Brunelle presented the paper, 

"Thermal BucJiling of Initially Stressed 
Thic.s. Plates". 

9/22-25 Euromech 147, Symposium on Damage Mechanics (Prof. 
Kri.mpl), Cachiin, France 

Professor Krempl presented the paper, "Life 
Prediction for Creep-Fatigue". 


TABLE 1-3 


COMPOSITES -RELATED f'tEETINGS/TALKS HELD AT RPI 
(April 30, 1931 through September 30, 1981) 


Topic Date 

The Use of Affine Trans- 5/17/81 
formations in the Analysis 
of Stability and Vibrations 
of Orthotropic Plates 


Speaker (s) 

Gabriel A. Oyibo 
Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
Engineering and Mechanics 
Rensselaer Polytechnic In- 
stitute 


Dornier and Composite Ap- 7/17/81 
plications 


Cy 1 indr ical ly Orthotropic 7/22/81 
Circular Plates with a 
Concentric Mole and with 
Properties Varying Radi- 
ally 


H. Gunter Helwig 
Design Specialist for Com- 
posites Materials and 
Structures Testing 
Dornier System Gmb II 

Christoph Muser 
Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
Engineering and Mechanics 
Rensselaer Polytechnic In- 
stitute 


Nonlinear Thermoolastic 
Effects in Surface Me- 
chanics 


8/6/81 Joseph M. Pitkin 

Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
' Engineering and Mechanics 
Rensselaer Polytechnic In- 
stitute 


Loaded Hole Analysis in 
Composite Materials 


8/13/81 Kenneth Bubeck 

Grumman Aerospace Corp. 


Finite Strain Transient 
Response of Thin Struc- 
tures 


9/10/81 Jose J. A. Rodal 
Kaman Avi-Dyne 


TABLE 1-4 

COMi’OSITES-RELATED VISITS TO RELEVANT ORGANIZATIONS 


by RPI Faculty/Staff/Students 
(April 30, 1981 through September 30, 1981 


Visited 


Date 


Purpose 


Naval Research 

5/4-5 

Prof. H. 

Lectured on acoustic 

Laboratory, 
Washington, DC 


Scarton 

emission techniques 
and demonstrated Sig- 
nal Processing of 
Acoustic Emissions 
from Composite Mate- 
rials 

Storl i nq-Winthrop, 
Rensselaer, NY 

5/12 

Prof. R. J. 
Dief endor f 

Gave lecture: "Anal- 

ysis of Experimental 
Data" 

Krautkramer- 
Branson, Inc., 
Stratford, CT 

7/19-24 

Prof. H. 
Scarton 

To discuss acoustic 
emission and ultra- 
sonic and NDE tech- 
niques and instruments 

IBM Research 
Laboratory, San 
Jose, CA 

8/4 

Prof. S. S. 
Sternstein 

Gave lecture: "Visco- 

elastic Properties of 
Composites" 

DepartJi\ent of 
Energy , 

Washington, DC 

8/23-25 

Prof. U. 
Scarton 

Took part in Review 
Panel activities on 
composite fly wheel 
program 

Battel le Columbus 
Laborator ios , 
Columbus, OH 

9/14 

Prof. M. S. 
Shephard 

Presented talk: "Fi- 

nite Element Model 
Generation" 
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Analysis”. This is the second year they were offered. Each 
course lasted a week, and they were run serially so that a 
participant might come for both on the seimc trip. The level 
of the material was planned to be particularly useful to man- 
agers of engineering structures activities who are involved 
in technical work but who may not have taken courses for 
several years. Because of the wide variety of special courses 
available throughout the United States dealing with this sub- 
iect matter, only rather unique aspects justify additional 
offerings in these areas. Our programs were planned to be 
unusual in the respect that "hands-on" experiences were in- 
herent and required for completion of each course. Use of 
hand-held programmable calculators, application of optimiza- 
tion programs using computer graphics, lay-up and cure of a 
simple part in graphite epoxy and testing tensile, shear and 
compression coupons to failure were all part of the first 
course. As part of the second course, workshop sessions were 
held in RPI's Computer Graphics Center. In these sessions 
the enrolled engineers and managers defined, "constructed" 
and analyzed their own finite element models on refresh gra- 
phics terminals. 

Faculty-staff and attendees for the courses were as 
follows: 

Composite Materials and Structures 

Program Director: 

R. G. Loewy 

Institute Professor 

Rensselaer Polytechnic Institute 
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Faculty: 

K. Judd Dicfcndorf 
Professor, Materials Eng. 
Rensselaer Polytechnic Inst. 

H. Gunther Helwig 
Design Specialist, Compos- 
ites Group 

Dornier System Gmb H 

S. Leigh Phoenix 
Associate Professor, Me- 
chanical and Aeronautical 
Engineering 
Cornell University 

Students : 

William Appleman 
Mechanical Engineer 
David Taylor Naval Ship 
R & D Center 

John T. Blaha 

Deputy Director of Flight 

Testing 

Edwards Air Force Base 

Huey D. Carden 
Aerospace Technologist 
NASA Langley Research Ctr. 

Dan Christmore 
Eng ineor 
Gates Lear Jet 

James P. Cycon 

Rotor Systems Design Analyst 

Skikorsky Aircraft 

John Emerson 

Structural Engineer Sales 
Manager 

Information Systems Den’ in 

Philip J. Granum 
Mechanical Engineer 
David Taylor Naval Ship 
R & D Center 

Thomas F. Grapes 
Senior Engineer 
Westinghouse Electric Corp. 


Stephen W. Tsai 
Chief, Mechanics and Sur- 
face Interactions Branch 
Air Force Materials Lab. 
Wright-Patterson AFB 

Dick J. Wilkins 
Sr. Engineering Specialist 
General Dynamics, Fort 
Worth 


David Hazel ton 
Development Engineer 
Oak-Fothergill , Inc. 

George J. Hewitt 
Structures Supervisor 
Heath Tecna Precision 
Structures, Inc. 

Alan Kemp 

Engineer Specialist 
Information Systems Design 

Patrick E. Kune 
Engineer 

Baltimore Division 

Martin Marietta Corporation 

Vladimir Roth 

Research Assistant Professor 
Rensselaer Polytechnic Inst. 

Charles R. Swats 
Aerospace Engineer 
Applied Technology Lab. 
Unites States Army Research 
and Technology Laboratories 

James Y, A. Wang 
Baltimore Division 
Martin Marietta Corporation 
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The Finite Element Approach in Structural Analysia 


Program Director: 

Mark S. Shephard 
Assistant Professor 
Civil Engineering 
Rensselaer Polytechnic Inst. 

Faculty: 

Richard H. Gallagher 

Dean. College of Engineering 

University of Arizona 

Mark S. Shephard 
Assistant Professor 
Civil Engineering 
Rensselaer Polytechnic Inst. 

Students : 

Gaston Beaulieu 
Structural Scientist 
Hydro-Quebec 
Varennes, Canada 

James Finkel 
Development Engineer 
General Tire and Rubber Co. 

Samuel M. Goldfarb 
Mechanical Engineer 
IBM, Poughkeepsie 

Aaron D. Gupta 
Mechanical Engineer 
United States Army Ballistic 
Research Laboratory 
Aberdeen Proving Ground 

Matt Harris 
Engineer 

Eastman Kodak Company 


Zbiegniew Karaszewski 
Mechanical Engineer 
Military Sealift Command, 
Atlantic 

Timothy Peters 
Project Engineer 
R & D Division 
Chamberlain Manufacturing 
Corporation 

Sanjay Shinde 
Senior Research Engineer 
Esso Research of Canada 
Calgary, Canada 

Tommy N. White 
Aerospace Engineer 
Naval Air Rework Facility 


Reaction to these offerings on the part of the partici- 
pants was again highly favorable, and we plan to make these 
offerings again in 1982. 
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4 . Sunmary 

In short « the NASA/AFOSR Composites Aircraft Program Is 
a multi-faceted program planned and managed so that scien- 
tists a^id engineers In a number of pertinent disciplines will 
Interact to achieve Its goals. Research In the basic compo- 
sition, characteristics and processing science of composite 
materials and their constituents Is balanced against the 
mechanics, conceptual design, fabrication and testing of gen- 
eric structural elements typical of aerospace vehicles so as 
to encourage the discovery of unusual solutions to present and 
future problems. In the following sections, more detailed 
descriptions of the progress achieved in the various compo- 
nent parts of this comprehensive program are presented. 
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II-A 


PART II 

CONSTITUENT MATERIALS 


A STEP IN DETERMINING TRANSVERSE FIBER PROPERTIES: 
CHARACTERIZATION OF HOMOGENEITY IN COMPOSITE MATERIALS 


.’RECEDING PAGE BLANK NOT FU MED 
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II-A A STEP IN DETERMINING TRANSVERSE FIBER PROPERTIES: 

C HARACTERIZATION OF HOMOGENEITY IN COMPOSITE MATERIALS 


Senior Investigator: 


R. J, 


Diefendorf 




Introduction 

Transverse properties are important characteristics of 
fibers to be used as constituent materials in structural com- 
posites. The small diameters of fibers ("Sm), however, makes 

it difficult to evaluate transverse fiber properties. This 
is especially true for the transverse thermal expansion of a j 

single fiber which requires accuracy oi: the order of several 
angstroms. An alternate approach is to evaluate various com- 

i 

posite laminates with differing volume fractions and plot the 
property versus fiber volume percent and extrapolate to 100% 
fiber volume percent. The implicit assumption is that it is 
far easier to measure volume fraction and dimensional change 
on a larger object and to accomplish the task with a greater 
certainty than for a single fiber. Therefore, the task trans- . | 

J 

forms into producing a reliable estimate for volume fraction ,1 

and a theory for transverse properties with anisotropic con- | 

stituents. As is well-known, the transverse properties of a j 

fibrous composite, unlike longitudinal properties, depend upon I 

J 

the fiber packing geometry. A subtle concept that has, as of | 

■a 

* . I 

yet, not been fully explored , is the effect of localized 

packing geometries on the properties. For example, a composite | 


A 

Holliday, Leslie (Ed) , Composite Materials , p. 11, Elsevier p 
Publishing Co., Amsterdam (1966) presents the only mention of | 
this phenomenon known to exist. 


rrtx jimmm ^ 
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that has an overall volume fraction of 65% may have an area 
with localized highs and lows of 72% and 50% respectively, 
while a second area (also with an overall volume fraction of 
65%) could have relative high and low volume fractions of 67% 
and 60% (Figure II-A-1) . The issue presented is to what de- 
gree, if any, does chis inhomogeneity affect composite proper- 
ties and performance. An extension of this concept is how can 
this index of homogeneity be quantified for use as a compara- 
tive measure and, if possible, to relate the index to the 
somewhat abstract (and often inadequately defined) concepts of 
poor and acceptable laminates. 

2. Status 

In order to obtain an estimate for the overall volume 
fraction of the fabricated laminates, the method of nitric 
digestion was employed. The samples were weighed, both in 
air and submerged in water, to determine specific gravity. 
Samples were then boiled in 100 ml of UNO^ (nitric acid) at a 
temperature of 177"C for 90 minutes and then drained and fil- 
tered through a porous ceramic crucible with acetone and 
water. The residual fibers were heated at 120“C for two 
hours and then cooled in a desicator. Fibers were then weighed 
and Boeing method MS 8-212A was used to calculate the volume 
fraction. The measurements are presented in Table II-A-1 for 
several of the many specimens tested. 
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3, Progress D trinq Report Period 
It has already been established that the accurate deter- 
ir.ination of volume fraction is critical to the work. It be- 
comes ajjparent that standard techniques alone, in this case 
nitric digestion, are not sufficient for the required accur- 
acy. It is cor this reason that two other methodologies are 
being evaluated to crosscheck the volume fraction, namely 
ASTM's photomicrograph grid technique and image analysis. To 
insure that only minimal discrepancies arise between the var- 
ious techniques, the latter two approaches utilize the same 
photomicrograph . 

A specimen. Number 18, was chosen from those listed in 
Table II-A-1 for further study. The sample was mounted in 
epoxy, polished and photomicrographs taken. In order to de- 
termine the volume fraction from the pictures via ASTM's grid 
techniq_e, a grid network consisting of lines spaced a quarter 
of an inch apart in both the horizontal and vertical directions 
on transparent mylar film was constructed. The grid area was 

TABLE II-A-1 

SPECIFIC GRAVITY AND VOLUME FRACTIONS 



USING 

NITRIC ACID 

DIGESTION 


Plate 

Number 

Specific 

Gravity 

Fiber 
Volume * 
Fraction 

Resin ^ 
Fraction 

Voids 

18 

1.6482 

68.3 

30.6 

1.1 

20 

1.5403 

68.1 

28.9 

3.0 

30 

1.4320 

49.7 

46.7 

3.6 

31 

1.5313 

69.4 

26.2 

4.4 


★ 

All values are percentages. 



Schematic Diagram of the Variatio 
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constructed so as to contain ten intersections on a side (10 
X 10) or 100 to the sheet. The grid was then laid at random 
on the photomicrograph and the number of intersections that 
fell on top of a fiber were counted, thus the total yields a 
sample value for volume fraction. This procedure was repeated 
ten times for each picture with the results averaged for an 
overall mear^ure of volume fraction. A listing of the data is 
given in Table II-A-2. It is conceivable that the grid size 

TABLE I I -A- 2 

COMPARISON OF VOLUME FRACTION DETERMINATIONS USING 
GRID AND IMAGE ANALYSIS TECHNIQUES 

Picture Number 



18-1 

18-2 

18-3 

18-4 

18-5 

18-6 

18-7 

18-8 

Image 

Analysis 

66.1 

70.0 

66.0 

63.9 

67.1 

67.2 

76.1 

71.7 

10 X 10 Grid 

68.9 

70.3 

65.9 

65.9 

65.5 

66.1 

72.5 

70.6 

5x4 Grid 

68.5 

66.3 

62.3 

67.5 

67.0 

63.5 

70.8 

69.3 


OVERALL STATISTICS 


Number 

Mean 

Standard Deviation 
Variance 


Image 

Analysis 

24 

68.4 

3.69 

13.03 


10 X 10 5x4 

Grid Grid 

80 80 

68.2 66.9 

4.39 8.46 

19.06 70.70 
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could affect the results of the analysis significantly. A 
5x4 grid was constructed and the analysis reperformed with 
the results being not too different from the original case. 

Image analysis was performed on a Bausch and Lomb FAS/I I 
image analysis system that utilizes a Data General micropro- 
cessor (Figure II-A-2) . Each of the eight pictures was re- 
solved for volume fraction with the results listed in Table 
II-A-3. It is important to point out that the values listed 
for each picture used in image analysis are the averages of 
three values, while for the case of the grids the values pre- 
sented are the averages of ten readings. An average overall 
reading is also presented in the table as is data for standard 
deviation and variance. The variance is of particular inter- 
est since it may be helpful in formulating a homogeneity index. 
The next phase was the subdivision of the main field into 
fields of equal area. The subfields were progressively divided 
into halves down to one-thirty-second of the original area. 
Volume fraction data was collected for each of the full, half, 
quarter, eighth, sixteenth and one-thirty-second fields. 
Averaged data for groups of ten measurements is presented in 
Table II-A-4 and Figure II-A-3. 

Each of the values of volume fraction in Table II-A-4 is 
an average of ten values. The values presented for standard 
deviation and variance result from the statistical analysis of 
these sets of ten measurements. The values of variance can be 
subdivided as follows, assuming that the contributing errors 




•‘S. IS* 






Figure II-A-2. 

Bausch & Lomb FAS/TI - Image Analysis System 





'OLCME FRACTION P'ES’JLTS FROM PHOTOMI . ROORAPH GRIO ANALYSIS 
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TABLE I I -A- 4 

VOLUME FRACTION DATA FROM SUBFIELD ANALYSIS 


32 

nds 

16'^'’“ 

jtho 

^ths 

Halves 

62.8 

52.7 

62.7 

63.5 

67.4 

68.6 

58.5 

58.1 

67.7 

65.4 

64.9 

64.9 

71.1 

66.9 

74.2 

63.0 

64.8 


66.0 

63.1 

63.0 

64.9 

68.3 


63.5 

69.3 

60.9 

62.9 



64.6 

65.5 

62.6 

66.4 



70.6 

65.4 

60.8 

60.8 



70.2 

71.2 

67.8 

65.6 



69.2 

76.5 

58.7 




53,2 

66.6 

60.1 




50.3 

70.1 

62.8 




64.2 

63.5 

65.4 




59.9 

60.1 

59.0 




70.6 

63.7 

64.5 




68.5 

50.5 

57.6 




66.2 

62.5 

74.3 





Full 

66.0 


Mean 

64.2 

63.9 

64.1 

66.4 

66.8 

66.0 

Standard 

Deviation 

6.325 

5.003 

1.842 

1.771 

2.616 

— 

Variance 

38.754 

23.465 

2.97 0 

2.353 

3,423 



VARIANCE 



VARIANCE VERSUS (AREA)"' 
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are independent. 


a 


2 

m 


”o * 


p/n 


( 1 ) 


Where: 

0 ^ is the sum of variances due to 
measurement, 

0 * is the variance due to random 
° events, 

0 ^ y is the variance due to the pixel 
P' ” count of the image analysis sys- 
tem divided by the number of 
fields. 


Due to the , term, is a function of the area of 
p/n m 

measurement and increases as the number of subfields increase. 
The data is plotted in Figure II-A-3 along with the "theoret- 
ical" curve determined by a manual, least squares fit assuming 
a slope of unity (the pixel errors are directly proportional 
to pixel count). The intercept of 0.15 implies that if a 
field of infinite dimensions were analyzed an infinite number 
of times, a variance of 0.15 would result. This variance is 

the 0 ^ term and is free from field size effects. Since a* and 
o o 

are known, can be determined. The importance of de- 
termining 0 ^ becomes apparent when a plot of volume fraction 
versus number of fields is constructed as in Figure II-A-4. 

In order to estimate the true sample-to-sample variance, it 
is necessary to correct for the contribution due to measure- 
ment error. This can be accomplished by means of the follow- 
ing equation which, once again, assumes that the errors are 
independent. 


NUMBER OF FIELDS 
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o 


2 



m 


( 2 ) 


Where : 

Is the overall variance for every 
measurement performed for a given 
subfield (i.e., for all 320 meas- 
urements in the case of 37nds, '60 
for 16ths, etc.)# 

Og is the variance due to the samples 
(unknown, to be solved for) and 

0 * is the measurement variance as de- 
fined in Equation (1) . 

After 0 ^ has been solved for, the standard deviation is 
s 

calculated by taking the square root. The "error” bars around 
the volume fraction data in Figure II-A-4 are one standard de- 
viation above and below the mean value. The increase in stan- 
dard deviation for small sample areas indicates an increase 
in the non-homogeneity of composite samples as the sample area 
is decreased. Comparison of this change in standard deviation 
with that obtained from the theoretical value associated with 
a uniform hexagonal or square array composite can be used as 
an index of homogeneity. 


4. Plans for Upcoming Period 
Future work will evaluate the importance of considering 
an index of homogeneity ir composite materials. Studies will 
continue on transverse properties in composites. Comparison 
of properties with the index of homogeneity will be made. 
Finally, transverse properties of anisotropic fibers such as 
carbons and aramids will be derived from composite data. 
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COMPOSITE MATERIALS 
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III-B FATIGUE IN COMPOSITE MATERIALS 
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PROPERTIES OF LAMINATES 
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III~A ADVANCED ANALYSIS METHODS 

Senior Investigator: E, J. Brunelle 

1. Introduction 

Because many parameters exist in a given composite (or 
orthotropic) plate problem, the established methiods of analy- 
sis are of little help in trying to answer the basic question 
"What are the general solution trends when I increment these 
variables and hold those variables constant?". A much clearer 
picture of how an individual orthotropic lamina (a plate) 
behaves has been made possible by previous wor)c under joint 
NASA/AFOSR sponsorship. This worlc has continued during the 
current reporting period by enlarging the scope of some pre- 
vious results and by introducing some new ideas. 

2. Status 

The prxiiCipal results of the last period were the follow- 
ing derivations, and some solutions, to a wide range of prob- 
lems associated with specially (and generally) orthotropic 
plates . 

a) The simplest inicromechanics description of 

D* = '^3s formulated with 

the main intent of showing 0 £ D* £ 1; this 
indeed was the case and additionally, light 
was shed on low fiber volume behavior - 

b) The generally orthotropic plate problem (and 
hence the generalized plane stress problem) 
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was transformed into rotated (thru an angle 
0) and affinely stretched coordinates such 
that only three coefficients remain in the 
formulation (L, D, R) and each of these co- 
efficients is only a function of D*, ^n /^22 
and the rotation angle G. 

c) The von K^rmdn-Rostovstev nonlinear plate 
equations were reduced to equations with 
only two parameters, D* and H*, where it 
was possible to show additionally that 

H* = H* (D*, f ) . 

d) Similarity laws were presented for a class 
of buckling problems which demand only that 
given number pairs must be invariants when 
proceeding from solution 1 to solution 2. 


3. Progress During Report Period 

Effort in a number of diverse areas has produced the 
following results (a through d) . 

a_) Attempts to find general solutions for the homogeneous 
problem associated with the general orthotropic plate have 
had very limited success; the literature is replete with in- 
teresting but specialized results. The first major stumbling 
block is that the equation to be solved, 






D ---- + D 




— - * °2, 
2^„2 26 


9‘^w 


+ D, 


9**w 




3y3x 


3 


0 , 


has six elastic constants whose ranges are not well-defined . 


If considered to be a single constant and if 


one divides them by any convenient constant (usually or 

** 

^ 22 ^* then four hybrid constants still remain. An appropri- 
ate affine transformation, reported on in the last progress 
report, reduces the problem to 




. 2D 


3x^ 9y* 
o -^o 




•*0 o 




where 

L, D ani R depend only on two parameters 

= 0 ^^*^0 = 0 rotation 

angle 0 . 

Xx wy 

Assuming complex exponential solutions w = e °e ° one finds 
p- + Ry^ + 2Dy^ + Ly + 1 = 0 

where 

y = — and where the "over-bars" will be 

10 

dropped in the sequel. 

Since complex roots are expected also for y, the factored 
form is left in terms of two quadratic expressions 
(y^ + a^y + b^)(y^ + a^y + b2) * 0 
where the and b^ (ot = 1,2) are given by the well-known re- 
sults 


However, their dependence on the rotation angle 0 from 
principal directions is known explicitly. 

* it 

They are called constants in most of the literature, but 
they are obviously parameters depending on the rotation 
angle 0 and their values at some reference 0 , usually 
0 = 0 . 
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bi + 

^ 1^2 ^^2 




R 

D 

L 

1 


Dividing the quartic by p** one finds 




2D — + 



0 


This is not usually of much help, but our constants L and R 
are related by 

L(0) = R(90 - 0) 

R(0) = 1,(90 - 0) 

D(0) = D(90 - 0) 

hence the second quartic may be written as (switching L and R) , 


+ R 


+ 2D 


L - + 1 
M 


0 


Therefore, if p^(i = 1,2, 3, 4) are the roots of the first 

quartic, then — are the roots of the second (above) quartic. 

^ i 

Furthermore, more information can be extracted from this pro- 
cess, and it is presented below as a new theorem of abstract 
algebra . 

Given the quartic in the two equivalent forms shown 


p” + Rp’ + 2Dp^ + Lp + 1 = 0 

and 

(p^ + a^p + b^) (p^ + a^p + b^) = 0 






then 


b^(R,D,L) = b^(L,D,R) 

with the corollary results 

b2(R,D,L) = b^CL^D^R) 

and 

aj^a2(R,D,L) = a^a2(L,D,R). 

Since bj^ and b^ will be positive for complex roots of the 
M -equation, it is seen that b^^ (or b^) has the range 0 to 1. 
Thus it is very quick and easy to find the a^ and b^ coeffi- 
cients. Furthermore, because of the L(6) = R(90 - 6) etc. 
relations, thr range of 9 to be examined is cut in half, 
which represents anot her large reduction in the operations 
to be carried out. 

This work directly impacts on all the literature con- 
cerned with finding homogeneous solutions for plate problems 
and all plane-stress/plane-strain problems. It currently 
appears that all of the useful but very clumsy procedures 
for solving complete plate problems (reorthogonalization ex- 
pansions, bi-orthogonal expansions of related first-order 
systems of equations, generalizations of the Fadle-Papkovitch 
techniques, etc.) can be expressed with much more clarity 
due to the above advances. 

b) Exact lower solutions bounds for the specially ortho- 
tropic plate can be obtained by letting D* = 0. Thus if we 
consider the vibration problem, the relation 
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^ 2D* 


+ 



Mw ■ 0 


( 1 ) 


+ + Mw - 0 

3x1 ay" 


( 2 ) 


becomes 

a'lw 

{** 

O •'O 

The first reaction is that one can hardly expect meaningful 

solutions for the D* > 0 case. However, this was exactly 

the case for the buckling solutions reported on previously; 

the D* " 0 solutions of versus a /b produced very useful 

o o o 

lower bounds, just as D* => i produced very meaningful upper 
bounds (the quasi-isotropic solution) . Those D* * 0 solu- 
tions were obtained by the limiting solution of the first 
equation (1) and not by the direct solution of the second 
equation (2) . Introducing 

in Equation (2) for a direct solution, one obtains 

^ ylVj^ _ „^2 xy = 0 (3) 


It is seen immediately that this is a separable equation, 
hence, exact solutions can be found for arbitrary boundary 
conditions on all four edges. Thus, for example, an exact 
solution can be obtained for the plate clamped on all four 
edges when D* = 0. Not only can the eigenvalue problems be 
solved (vibration and buckling) , but the static deflection 
problems also can be solved by expansions in terms of the 
eigenfunctions (mode shapes) of the associated homogeneous 
solution (this is a standard technique of mathematical 
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physics) . Furthermore I the dynamic response problems for 
D* = 0 can also be obtained in terms of generalized coordi- 
nate** weigh ted eigenfunction expansions (again a standard 
technique) . Separating Equation (3) one has 


-f )«'* - Mu)** X = 


0 and - )c**Y 


where the separation constant, for illustrative purposes, 
has been chosen as -k** . It is noted that the separation 
constant can be +k‘*, -k** or 0. The correct choice is that 
value which allows the boundary conditions of the problem to 
be satisfied. 

As an example of this procedure, it can be verified that 
the exact solution for the non-dimensional frequency of the 
SS-C-SS-C plate with D* = 0 is 


2 ^ = 
mn 


mb 

o 

4 

r 

2n - 1 

a 

f 

4 

o J 




where 


mn mn( 


Furthermore, for small values of 2D* = e one may write 


8 w , _ 9^w . 9^w , ^ 

+ £ + = q(x^,y^) 


9x‘ 


9x"9y^ 9y‘ 

o -^o 


and using regular perturbation expansion techniques, one may 
assume 


w = E w (x ,y ) 

n o'^o 

n=o 
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so that the sequential set of equations 






3x‘ 


ay** 

•^o 




a ^w. 


a^w. 


a-w. 




ay“ 

•'o 


o ^ o 


a**w, 


k+l 


a^w. 


k+l 


a^w. 


ax‘ 


3y. 


ax^ay^ 
O o 


o 1 

may be solved to any given k value. Thus the stage is set 
for solving not only static problems (the above development) 
but also stability and dynamics problems for small e. It is 
noted that one is again able to specify arbitrary boundary 
conditions for all these problems. 

c) D* is the fundamental plate parameter for the specially 

orthotropic plate. Its behavior, via simple micromechanics, 

has been displayed versus the fiber volume for constant 

values of E^/E , and v . From that perspective, much in- 
t m t m 

formation about the incremental behavior of D* can be obtain- 
ed. An even simpler way of looking at the incremental be- 
havior of D* is to express dD* in terms of the macroscopic 
definition of D*. Thus, since 


D* = 


“12 




it can be verified that 




= e 
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where the identity D* = gD* + (1 - e)D* has been Invoked. 
It is thus seen that for a given incremental value of the 


plate constants, the 


d066 


D 


66 


terra yields the largest contribu- 


dD* 

tion to the expression since its coefficient (1 - e) is 

dDi 2 

larger than the e coefficient of or larger than the -•» 

dD, , dD,n 12 

coefficient of both -k and . 


11 


22 
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With this established, -r-^ and 

^o 


may be found, and as 


an example, the results for a simply supported specially 

orthotropic plate are given below (where the similarity num- 
b m 

ber z = has been used) . 

“o 


dk 


dt^2 

mn 


D* 


D* + 1 
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z“ + 2D*z^ + 1 
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d z 


+ D 


.dD* 
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d) The similarity rule., for vibration of a specially ortho- 
tropic plate have been derived and are being used in a nearly 
completed Masters' Thesis to compute the frequency spectrum 
^ral ' ^ ~ 1,2,3,... versus 3^/1)^ for the D* range of 0 to 1. 

Six sets of boundary conditions were chosen for the y =0 

o 

and b^ sides, the x^ = 0 and a^ sides being simply supported. 

The use of the affine transformation technique and the 
simile rity rule have permitted, for the first time, enormous 
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eimounts of frequency information (including important trends) 
to be presented in an extremely compact fashion. Letting 


»o, 

mb. ' 
o 


the similarity rules can be written as: 




fD*l 




and 


mn 1 


Q 

2 

1 - 

fD*l 

2 

-1^ 



(♦2)=D* 


X 

.’“I; 


f 


IH 2 


when D* = 0. When D* = 0 it can be verified that a single 
rule exists, which is 


<“mn> I 


- 




-1 + 


mn 2 


(4'2)'* 


(4>i) 
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4. Plans for Upcoming Period 
Results are being compiled to furnish a general overview 
of the affine transformation process and the similarity laws 
for orthotropic structures and the frequency spectrum of 
specially orthotropic plates using the similarity laws for 
vibration- This latter work is that of graduate student 
Chikuang Chen. 

At this time the new theories presented have not had 
enough supporting problems worked out to properly display 
their usefulness. As seems needed, specific problems v’ill 
be worked to illustrate points of interest. 
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A recent composite box-beam wing theory due to Weisehaar 
will be re-cast in its affine plane form to provide results 
of greater generality and will be extended to include simple 
non-symmetric flight maneuvers that provide large bending- 
twist-extension effects which yield opposite trends on the 
left and right wing. A senior student (Mr. Don Argintar) 
will be doing a senior project next semester based on this 
rederivation and extension. Forward sweep, aft sweep and 
oblique sweep configurations will be considered. 


5. Current Publications or Presentations by 
Processor Brunelle on this Subject 

"Thermal Buckling of Initially Stressed Thick Plates" 

Presented at the 4th Failure Prevention and Relia- 
bility Conference (Joint ASME Vibrations, Design 
Automation, F-P, R Conference) , Hartford, Connec- 
ticut, September 20-23, 1981. 


Weisshaar, Terrence A., "Aeroelaatic Stability and Perform- 
ance Characteristics of Aircraft with Advanced Composite 
Sweptforward Wing Structures", AFFDL-TR-78-116 , Sept. 1978. 
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III-B FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 

1. Introduction 

The deformation and failure behavior of graphite-epoxy 
tubes under biaxial (tension* torsion) loading is investi- 
gated as part of the research of Ph.D. candidate T.-M. Niu. 
The aim of this research is to provide basic understanding 
and design information on the biaxial response of graphite/ 
epoxy composites. 


2. Status 

Static tests on thin-walled graphite/epoxy (±45)^ tubes 
showed that the axial tensile and compressive strength ate 
approximately equal. However, the ultimate strength in tor- 
sion depends strongly on the direction of twist. Axial fa- 
tigue tests under completely reversed, load-controlled con- 
ditions showed a comparatively poor fatigue performance. 

Also, there was an influence of frequency on fatigue life. 

3. Progress During Report Period 

The major emphasis was on preparatory work for fatigue 
tests under combined in-phase and out-of-phase loading. A 
total of ninety tubular specimens were made and machined. 

They are ready for testing. 

During this reporting period, the project used a new 
463 Data and Control Unit using a PDP 11/23 computer and a Tek- 
tronix 4025 terminal in conjunction with the MTS tension-torsion 
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servohydraulic system. Driver programs were written and de*> 
bugged for biaxial in-phase and out-of-phase loading. Diffi- 
culties were encountered initially with the MTS software (RPI 
received the first 463 Data and Control Unit) which are now 
resolved. The programs have been tested with and without 
graphite/epoxy tubes. The driver programs and the data ac- 
quisition system are now fully operational. 

A technical report entitled "Graphite/Epoxy (^45)^ 

Tubes: Their Static Axial and Shear Properties and Their 

Fatigue Behavior under Completely Reversed, Load -Control led 
Loading" was written and is the first item in Section 5 be- 
low. It has been sutanitted for publication. 

4 . Plans for Upcoming Period 
Fatigue tests under completely reversed, load -control led 
conditions will be performed for biaxial in-phase and out-of- 
phase cycling. 


5. Current Publications or Presentations by 
Professor Krempl on this Subject 

1) "Graphite/Epoxy (±45)g Tubes: Their Static Axial and 

Shear Properties and Their Fatigue Behavior under Com- 
pletely Reversed Load Controlled Loading" 

Published as RPI Report No. RPI CS81-2, September 1981, 
(with T.-M. Niu) . 

2) "Preliminary Experience with the Instron Biaxial Extenso- 
meter" 

Presented at the ASTM E9.08 Committee Meeting, Phoenix, 
Arizona, May 13, 1981. 

3) "Life Prediction for Creep-Fatigue" 

Presented at the Euromech 147, Symposium on D<mage 
Mechanics, Cachan, Prance, September 1981. 
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III-C MOISTURE AND TEMPERATURE EFFECTS ON THE MECHANICAL 
PROPERTIE3 or LAMINATES 

Senior Investigator: S. Sternstein 

1. Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of the matrix resin. To date, the re- 
search has involved the precise viscoelastic characterization 
of epoxy neat resins, interlaminar failure of composites and 
the inhomogeneous swelling of and the effects of moisture on 
composites. 


2. Status 

In the previous report, the dynamic mechanical proper- 
ties of carbon/epoxy laminates were presented. It was shewn 
that the dispersion characteristics of neat resin and in situ 
resin were equivalent and that time-temperature superposition 
was obeyed in the glass transition region of the epoxy. Us- 
ing the centre-symmetric deformation (CSD) testing method 
described previously, we have investigated the effects of 

water on the dynamic mechanical properties, as reported here 

* 

for fiberite laminates . 


Fiberite prepreg HYE-1048AIE was used to prepare laminates. 
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3. Progress During Report Period 
Repetitive dynamic mechanical tests have been performed 
on dry laminate s'.mples, and it is found that the glass tran- 
sition temperature increases slightly from run to run until 
a final, fully cured Tg of 160 ‘'C is obtained. If a postcured 
laminate sample is subjected to moisture, then it is found 
that the glass transition temperature varies directly with 
moisture content as measured by weight gain. This is shown 
in Figures III-C-1 and 2 which give the in-phase and out-of- 
phase components or the dynamic modulus at 6.5 Hz, respec- 
tively. The glass transition is shown versus weight gain 
(water) in Figure III-C-3. The effects shown in Figures III- 
C-1 and 3 are reversible, in that a sample which is dried 
will recover its original glass transition temperature. 

Laminate samples which have not been previously post- 
cured display quite different effects when exposed to mois- 
ture. The initial exposure to moisture results in a decrease 
of glass transition temperature which is not recovered after 
drying. Typically the decrease in Tg is five to ten degrees, 
depending on the moisture treatment. Subsequent dynamic 
mechanical tests on wet samples indicate that moisture inter- 
acts strongly with the postcure Icinetics. We are currently 
investigating these kinetics by performing dynamic modulus 
tests versus time, at constant temperature. It appears that 
even after the sample has dried there is still a residual 
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effect on postcure kinetics. These results are being pre- 
pared for publication and will be available shortly. 

Finally, we have investigated the dynamic behavior of 
both wet d dry neat resin samples of Narmco 5208. The wet 
samples were conditioned at 62*C, 100% R.H. vapor for about 
2 years which was the time required to reach weight gain 
equilibrium. The in-phase and out-of -phase moduli! at 6.5 Hz 
are shown for the dry sample in Figure III-C-4 and the wet 
sample in Figure III-C-5. The broad out-of-phase modulus 
peak (M'') in Figure III-C-5, extending from 160®C to 240“C, 
should be compared with the peak at 240“C in Figure III-C-4. 
This broadness of the wet glass transition peak may be di- 
rect rheologocal evidence for the presence of structure in 
the resin. Pr-^viously, other investigators have suggested 
that an amorphous structure exists in thermosets, based on 
electron microscopy of fracture surfaces. Since Figures 
III-C-4 and 5 represent bulk behavior, not surface behavior, 
these results may be of considerable utility in probing the 
structural features of epoxy resins. 

4. Plans for Upcoming Period 

During the next period we plan to continue research to 
characterize the effect on dynamic mechanical properties of 
moisture interaction with neat resins and composites. Addi- 
tional studies on resin cure rheology and inter-laminar 
creep will also be initiated. 
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5. Current Publications or Presentations by 
Processor Sternstein on this i^ubject 

1) "Properties of ^ situ Epoxy Matrices" 

Presented Invited Lecture at the U. S. A. - Italy 
Joint Symposium on Composite Materials The Role 
of the Polymeric Matrix in Their Processing and 
Structural Propertiej, Capri« Italy, June 15-19, 

1981. 

2) "Viscoelastic Properties of Composites" 

Presented at the IBM Research Laboratory, San Jose, 
California, August 4, 1981. 

3) "Viscoelastic Characterization of Neat Resins and Compos- 
ites" 

Presented at the National A.C.S. Meeting, New York 
City, New York, August 27, 1981. 


58 


III-D NUMERICAL INVESTIGATION OF MOISTURE EFFECTS 

Senior Investigator: M. Shephard 

1, Introduction 

This analysis builds on the theoretical developments of 
Professor S. Sternstein (reported in progress reports dated 
December *79 and June '80) in which the fully-coupled, non- 
linear thermomechanical equations of inhomogeneous swelling 
of composites in the presence of moisture and temperature 
are analyzed. Here the one -dimensional case of a single 
fiber in an infinite matrix is being extended to two-dimen- 
sional, multiple fiber cases. 

2. Status 

As discussed in the previous progress report, moisture 
effects have been introduced into the problem through a non- 
linear constitutive relation. The principle of virtual dis- 
placements is employed to develop a nonlinear matrix equation 
using displacement-based finite elements. These nonlinear 
matrix equations are then solved using quasi-Newton methods. 

3. Progress During Report Period 

Effort has been continued during the last reporting 
period by graduate student Frida Lumban-Tobing in testing 
the program on some sample problems and establishing addi- 
tional features in the program. The first problem considered 
was that of a single fiber embedded in an epcxy matrix. This 
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problem was considered because « as has been noted, it is a 

one*dimensional problem that has been solved under the super- 

* 

vision of Professor S. Sternstein . The constitutive equa- 
tion employed in this one-dimensional analysis is written in 
the principal direction as 


2t 

o.V. = aUT — - 

1 ) T T 


1 kL 


3AV 

^ * In. 


12 


+ i”.(T 

ail ^^2 ^ 


- X^) 1 

^ 3^3 Xj 


( 1 ) 


where ; 


= principal stress 

I ,I,,I = invariants of the deformation tensor 

L a! I 

W A W(Ij ,I ,,I == the w^rk function 

X^ = 1 + 2e . ; c. . = the principal strain tensor 

1 3 1 f t- 

= partial molar volume- of the solvent in 
polymer 

= molar volume of pure solvent 

= volume change as a result of mixing poly- 
mer and solvent, per unit volume of poly- 
mer 

Tij = number of moles of solvent absorbed, per 
unit volume of polymer 

a = In (aj); ai = the activity of the solvent 
in the polymer 

R = gas constant 

T = absolute temperature 


Herbert, P. F., "Moisture Induced Swelling in Graphite/Epoxy 
Composites", Masters Thesis submitted to the Graduate Facul- 
ty, Rensselaer Polytechnic Institute, 1981. 
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After rearranging* this expreesion can be fit into the non- 
linear constitutive relationship employed in the program as: 


“ij ■ *«ij * * ‘^=ikSj 
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The basic finite element idealization employed consists 
of a fiber of radius r embedded in a 6r radius cylinder of 
matrix material as shown in Figure III-D-1. The first check 
analysis done was a free swelling analysis where both the 
outside and fiber boundaries are allowed to move. Figure 
III-D-2 shows the displaced mesh for this analysis. As would 
be expected, the stresses throughout the swollen medium were 
essentially zero. The second check analysis was for the case 
in which all the boundaries are fixed. In that case, the 
results were as expected, with the matrix in uniform compres- 
sion . 

The boundary conditions for the analysis of the reference 
on the proceeding page are defined as fiber boundary fixed 
and outer boundary free. This problem was solved for various 


Figure III-D-1. 


Finite Element Mesh Used in Single Fiber Analysis 



Figure III-D-2. 


Displaced Shape for the Free Swelling Case 
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moisture contents. Some numerical ill conditioning problems 
due to the log terms were experienced. These were allevi- 
ated, however, by adding features to the standard quasi- 
Newton procedures. The stability of the quasi-Newton approach 
was improved by adding line searches which improve the di- 
rection of the steps taken. In addition, when the size of 
the residual vectors exceeded a pre-set threshhold, one step 
of standard Newton iteration was performed. 

The solutions obtained compared very well with those 
presented in the reference on page 57. Figure III-D-3 shows 
the displaced shape for 80% humidity, while Figure III-D-4 
shows the principal stress contours at the same humidity. 

4. Plans for Upcoming Period 

Development of the nonlinear program for moisture ef- 
fects should be completed in the next reporting procedure 
for the multi-fiber, two-dimensional use. A series of analy- 
sis studies will be run. These studies will consider various 
composite configurations and conditions including fiber 
spacing, epoxy properties, moisture contents and loading con- 


ditions . 
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III-E NUMERICAL INVESTIGATION OF THE MICROMECHANICS O F 
COMPOSITE FRACTURE 

Senior Investigator; M. Shephard 
1. Introduction 

To understand the mechanisms of failure in composites 
it is necessary to develop insight into the micromechanical 
behavior, including interactions between matrix and fibers 
=^s the load is increased from zero co that corresponding to 
failure. Investigations of these phenomena, either experi- 
mental or nu.'ierical, are difficult. The purpose of this 
project, being carried out by graduate student Nabil Yehia, 
is to develop a nonlinear finite element analysis capability 
for performing numerical investigations of the micromechan- 
ical failure of composites. 


2. Status 

To gain some feel for the processes involved in compos- 
ite fracture at the micromechanical level, a step-by-step 
analysis of the failure of a simple three-fiber system was 
planned where one fiber is initially broken. Each step would 
consist of a linear analysis where the load would be increased 
until the next model modification would be called for. In 
addition to demonstrating how difficult it would be to per- 
form 2 step-by-step analysis, the results were expected to 
indicate that one must track the formation and growth of a 
crack, consider material nonlinearities, and in some very 
localized region, consider large strains. 
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Based on previous published works, it was anticipated 
that to be able to adequately model the desired behavior we 
would have to develop a nonlinear finite element program 
with the following features; 

(a) nonlinear constitutive relations, 

(b) nonlinear strain displacement relations and 

(c) crack modeling abilities. 

If only the first two forms of nonlinear behavior were re- 
quired it would be possible to employ available nonlinear 
finite element codes. However, crack modeling ability is 
not available in existing programs. Thus, the design and 
development of a special nonlinear code was indicated. 

3. Progress During Report Period 

In addition to a survey of the applicable literature, 
tht basic structure of the required FEM program has been de- 
signed and many of the procedures to be included selected. 
Coding of the first set of features has begun. The follow- 
ing paragraphs summarize approaches to be incorporated in 
the program. 

The matrix material tends to be nonlinear for most values 
of scrain, thus a constitutive relation based on a strain 
energy function was selected because of its ability to repre- 
sent general nonlinear relationships. For the cas 'f two- 
dimensional isotropic materials, for example, the relation- 
ship for the direct stress terms take the form (see, for 
example. Reference III-E-1) ; 
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and the shear inq stress is of the form: 

ii - R o + R e ’ 

The ineliTstic portion of the behavior will initially be 
lumdled employinq standard yield criteria (Ref. IIT-K-2) . 

Larqo strains will be handled usinq the standard non- 
linear strain-displacement relations (neforence in-C-3). 

The ability to model crack initiation and qrowth re- 
quires several features that must be built into the proqram. 
.'irst, one must be able to determine when a crack will initi- 
ate. Then one must determine when a crack will |)ropaqate, 
in whicli diii'ction it will qrow and how far it will qo. This 
in turn requires the ability to chanqc the qeometry of the 
finite element model. The methods to be employed to attack 
this problem will rely heavily on fracture mechanics concepts 
and the use of special sinqular finite elements. In particu- 
lar, the fracture criteria to be employed include: 

a) For cracks inside one material, the strain energy 
density factor S-criterion will be used (Ref. 

1 T I -E-4) . 

b) For cracks at a material interface, the S-criter- 
ion (Ref. III-E-4) and a maximum bond strength 
criterion will be used. 

c) In computing the stress intensity factors, K, 
the displacement method will be used (Ref. 111- 
E-S) . 
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d) The crack tip region will be modeled by a layer 
of six-noded singular elements (Refs. III-E-6 
and 7) followed by one or two layers of eight- 
noded transition elements (Ref. III-E-8) . 

e) The mesh of the region (s) in which the crack is 
embedded will be regenerated after each crack 
preparation. The scheme for doing that has not 
yet been worked out. 

4. Plans for Upcoming Period 
Efforts in the coming reporting period on the micro- 
mechanical failure of composites will be concentrated on im- 
plementation and initial testing of the various features 
being built into the analysis package. Attempts will also 
be made to devise a general procedure for handling mesh 
modifications as required by crack grow during load increases. 


5. References 
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PART IV 

GENERIC STRUCTURAL ELEMENTS 


IV-A 727 ELEVATOR ACTUATOR ATTACHMENT RIB 
IV-B LOCKHEED L-1011 ENGINE DRAG STRUT (CAPCOMP II) 
IV-C MECHANICAL JOINTS IN COMPOSITES 
IV-D COMPUTER SOFTWARE DEVELOPMENTS 
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INTRODUCTION TO PART IV 

The redesign of flight-critical structural components 
currently in service as metal parts on existing aircraft was 
included as part of RPI's composite materials and structures 
program to assure relevance and realism to the project's di- 
rection and to contribute to the nation's pool of innovative 
structural design concepts. Two such components were se- 
lected in cooperation with major commercial airframe manu- 
facturers, Boeing Commercial /Mrplane Co. and Lockheed Cali- 
fornia Comp ;iy, and with advice from NASA. The first, the 
7 27 aii. plane's elevator actuator attachment rib, has com- 
pleted the cycle planned at the project's inception. That 
is, several alternate redesign‘d were attempted? one was se- 
lected, analyzed and finalizes; and two specimens were fab- 
ricated and tested to ultimate load. The second, a redesign 
of the wing-mounted engine drag strut for the L-1011 airplane 
has been carried to the point where (a) the usefulness of 
certain basic research performed earlier at RPI under the 
NASA/AFOSR program was made very clear and was called on in 
the design process, and (b) several areas of research in 
which results are needed for such a redesign - but are not 
now available - have been identified. 

As a matter of fact, a number of basic research projects 
motivated by the early work on the 727 elevator rib, have 
been reported in these project reports for some time. Some 
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of these have been completed, others are reported here in 
connection with the conceptual design research on the two 
specific airframe components, as discussed above, and still 
others are planned for the future as having the potential 
to contribute significantly to classes of generic Composite 
Aerospace Structures, as the specific aircraft structural 
component redesign research is phased out. Progress on 
projects typical of the original program and projects of 
generic research, as will be pursued in the redirected pro- 
gram, are both reported in the sections that follow. 
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IV- A 727 ELEVATOR ACTUATOR ATTACHMENT RIB 

Senior Investigators i D. Goetschel 

R . Loewy 
H. Scarton 

Tnis program element involved a structures demonstra- 
tion project using the actuator attachment lugs and rib of 
the 727 elevator, paralleling that of NASA and its aerospace 
contractor, the Boeing Commercial Airplane Company. Our de- 
sign, fabrication and test effort emphasized designs using 
advanced composite construction for the purpose of minimiz- 
ing the weight of the structure, but on a scale consistent 
with the university context and funding level. The staff of 
RPI is very grateful to the Boeing Company and its engineers 
for their wholehearted support of this work at RPI. 

One of two different designs was reported previously 
(July and December, 1973) as suitable for replacing the 
largely metal attachment produced by Boeing. This so-called 
"Berg Design" makes use of quasi-isotropic, graphite/epoxy 
laminates and was selected for further analysis, redesign, 
fabrication and testing. The conventional aluminum alloy 
elevator of the Boeing 727, the new version of the elevator 
redesigned by Boeing in graphite epoxy, and Berg's design 
for the RPI graphite/epoxy Nomex honeycomb elevator rib are 
all described in earlier reports (July, 1981) . In each case 
the attachment was designed to carry loads up to 19,000 
pounds and for directions of load varying as the elevator 


74 


rotates over an angle of 28 degrees from the Cull-down to 
the full-up position. 

The first specimen of the Berg-designed, composite 727 
aircraft elevator actuator attachment rib was fabricated by 
undergraduate student Francis DeTaranto and was tested to 
failure on December 12, 1980. As reported in the semi-annual 
report dated July, 1981, this was a completely successful 
test in virtually all respects. The design failed at approx- 
imately 120 per cent of design ultimate load, failures oc- 
curred close to where they were expected and failures were 
multiple in nature. The major failure appeared to have 
initiated near the leading end of the tension flange (i.e., 
the flange connecting the rib to the elevators' upper-surface 
skin) . 


1. Construction Details of the Final 727 Elevator Rib 
During the current reporting period a second and final 
test specimen was fabricated by undergraduate student Frank 
DeTaranto and graduate student Shu-Ping Chen. Four adjust- 
ments were made in the detailed design of and processes used 
to fabricate the second part. The changes were intended to 
improve the part's behavior in the area where failures initi- 
ated in tests of the first rib. 

First, "wrinkles" were observed to have developed in 
fabrication along the inner radius of some of the flange-web 
junctions, running parallel to it. It appeared, initially, 
that the cause of the "wrinkling" was the bunching-up of 
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fibers from shrinkage during curing or from vacuum bag drag- 
ging outer layers into overlapping as pressure was applied 
and the bag tightened. None of the attempted remedies « how- 
ever, had a significant effect. 

Photomicrography finally revealed that the real cause 
of the "wrinkling" was the separation of zero degree fibers 
in the surface ply, and their piling-up on either side of a 
"chasm" to cause what only appeared to be a wrinkle. Figure 
IV-A-l-a is a typical micrograph of the "wrinkle" with 0®/ 
outer plies (i.e., parallel to the corner). Reordering of 
the stacking sequence cured the problem. In the second part, 
the surface ply was oriented at 90® (i.e., perpendicular to 
the axis of the corner) . This also increased the local 
bending strength and stiffness of the flange. 

To describe the second design change, one must call at- 
tention to the build-up of the end of the flange on the first 
actuator rib (Figure 1-20 in the July, 1981 report, repro- 
duced here as Figure IV-A-l-b) . This build-up increased the 
stiffness of the flange at the flange leading edge, allowing 
me e load to be picked up by the flange in the region of the 
first fastener. This is where failure occurred. Thus, on 
the final actuator attachment rib, this build-up was elimi- 
nated . 

The third modification was spreading Hy-lock fastener 
locations away from the web as the chordwise position of 
these connectors between skin and flange moves from trailing 


diai.ii< 


MiliilbllfiilllM 




m 


Fillet Crease 
flAGN! FI CATION SOX 





Failed Section and Strain Gage 


78 


to leading edge. The actuator rib web becomes thicker near 
the leading edge of the rib to carry the more concentrated 
loads associated with the actuator attachment lug. The 
fillet of the flange thus is dicker in those areas, result- 
ing in a sharper taper in flange thickness from fillet to 
edge. This, in turn, requires a wider fastener spacing to 
provide a flat "seating” surface for the fasteners. 

The fourth and final modification was the use of washers 
under the Hy-lock fasteners to improve the uniform distribu- 
tion of stress between fastener and composite. 

2. Test Procedure and Results for the Final 727 Elevator Rib 

Oi September 11, 1981, the revised Boeing 727 elevator 
actuator attachment was tested. Figure IV-A-2-a shows the 
schematic of the test apparatus, which had been modified 
from that used in testing the first specimen by including 
spherical-ball bearings at the pin connections (points A in 
Figure IV-A-2-a) and load-carrying transition pieces at the 
lower section of the loading fixture (point B in the figure) . 
Figure IV-A-2‘-b shows the elevator attachment placed in the 
Instron 1333 machine whose maximum capacity is 252 kN (56,500 
pounds) . A continuous monotonically increasing loading with 
constant load rate was used. In the final test run, load 
was supplied at 0.135 kN/s (30 pounds/sec) to the specimen. 

The instrumentation used in the test is shown in Figure 
IV-A-2-C. It consists of an FX Data Logger Strain Monitor 
with four 45® strain rosettes and seven unidirectional strain 
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gauges. A Physical Acoustics Corporation Model 3400 four- 
channel acoustic emission monitor supplied for the test 
through the courtesy of the Physical Acoustics Corporation 
was also used. RPI graduate student Chi-Min Chang and Mr. 
John Crowell from PAC supplied strain gauge and AE instru- 
mentation support, respectively, both during and after the 
test. A personnel protection barrier is also shown in Figure 
IV-A-2-C. Table IV-A-2 lists the acoustic sensors used in 
the Lest. 

In order to place strain gauges and AE sensors as close 
as possible to the ultimate failure location, a simplified 
two-dimensional finite element analysis was performed using 
a finite element mesh which represented thickness variations 
as shown in Figure IV-A-2-d. Figure IV-A-2-e presents the 
overall maximum principal stress contours which resulted 
from this analysis. The results of a pre-test study of the 
lug, shown in Figure IV-A-2-f for maximum principal stresses 
and shown in Figure IV-A-2-g for specimen distortion in the 
vicinity of the edge of the upper pin hole, suggested that 
this was likely to be the region of ultimate failure. On 
this basis the AE sensors and strain gauges were located as 
indicated in Figure IV-A-2-h. The AE sensors were calibrated 
using five clicks on a 0.5 mm 2 H hardness Pentel pencil and 
then by breaking the lead held at an angle of 45®. These 
im.pulses have consistent amplitudes of 8 0 dB. 

The actuator was then loaded to 2,000 pounds and all 
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Fxgure IV-A-2-d. 

Variable Thickness Homogeneous Finite Element Spe ;imen 
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Figure IV-A-2-e. 

Finite Element Maximum Principal Stress Distribution 
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Figure IV-A-2-f. 
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Figure IV-A-2-g. 

Finite Element Simulation of Final Strained Shape 
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strain gauges Interrogated. Although some asymmetry in the 
strains was reported, a decision was made to proceed with 
the test. 

The load was removed and the Instron was readied. An 
initial ramp loading of the specimen was carried out, to a 
load of 55.9 kN ( 1 2 ,600 pounds) (Figure IV-A-2-i) . Because it 
was discovered that a lower loading rate was needed, the 
test was halted at this point and the load removed. It was 
confirmed in poLt-test analysis of the AE data that, during 
the first loading, over ten large amplitude events (exceeding 
90 dB) were observed (Figure IV-A-2-j). The loading rate was 
reduced by a factor of 10 to 0.135 kN/s (30 pounds/sec). 

The test was then restarted with all AE data channels 
reinitialized to zero. The test continued until ultimate 
failure occurred at 721 seconds, at a load of 96.2 kN (21,630 
pounds) or 14% over the design load of 85.5 kN (19,000 
pounds) . 

During the test cumulative counts versus time were moni- 
tored in real time (see Figure IV-A-2-k) . A pronounced 
Kaiser effect can be seen in the second run, where no counts 
were observed until the load exceeded 53 kN (12,016 pounds, close 
t"> the maximum established during the first ramp loading. 
Alti'vough there are a number of possible points of critical 
activity wherein the slope of the curve of number of counts 
versus load increases precipitately, two points appear to 
have noteworthy critical activity; one at 66.56 kN (15,015 
pounds) and another at 84 kN (18,700 pounds). Post-test 
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Figure IV-A-2-i. Specimen Load Versus Stroke Curve 
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examination of the eunplitude distributions from th<=* second 
data run with a slower loading rate shows a broader distri- 
bution of events with far less fiber breakage (i.e., fewer 
large amplitude events} . The final amplitude distribution 
at failure (shown in Figure IV-A-2-il) indicates a curious 
build-up of counts at the 65 dw level. The 65 dB amplitude 
events (totaling 144) were observed throughout the test and 
may well be directly related to the specific failure mecha- 
nism. Also, there appears to be a reasonable correlation 
between the eleven sudden increases in AE counts (Figure IV- 
A-2-k) with the eleven load-dvop points that are seen (Figure 
IV-A-2-i) in the discontinuous slope region of the load- 
stroke curve. 

Examination of the strain gauge data shows a lack of 
symmetry in the readings between the two sides. At four of 
five comparison stations, side two had readings about 50% 
higher than those on side one. The fifth station showed side 
one readings 300% higher than side two. Examination of the 
two sides of the test rib shows significant differences in 
the chara -eristics of this fifth position, from side one to 
side two. These differences apparently came about from the 
combination of individually small inaccuracies in strain 
gauge placement and small variations in laminate stacking 
boundary positions between the two sides, in a region of 
rapid cross-section variations. The accumulation of these 
effects are believed to be the most likely explanation for 
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this 300% disagreement. The imbalance at the other four 
gauges is probably due either to additional asymmetries in 
the part itself or an imbalance in the application of the 
loading. There is little consensus on which is the correct 
explanation. Contact with Boeing acquired the information 
that ir. aircraft practice, care is taken to position loading 
pins ccwtrally between lugs. Spherical bearings are used to 
prevent moments from iuadvertently entering. In spite of 
this, safety margins of 20% or more are added to account for 
asymmetries. The failure actually occurred in our test on 
side one, which had the generally smaller measured strains 
(Figure IV-A-2-m) . For example. Figure IV-A-2-m shows the 
u.iidirectional strain values for vertical strain at points 
E and EP (Figure IV-A-2-h) . Figure IV-A-2-n shows the prin- 
cipal strain values and other selected strains at 707 seconds 
into the test. The load at this time was 94.3 kN (21,210 
pounds). Using the material properties E = 40 GPa (5.83 Mpsi) 
and V = 0.28, the principal stress values were computed from 
the principal strains. (Note that this condition is very 
near the ultimate failure that occurred only 13 seconds later.) 
Marked asymmetries were observed to exist, yet the observed 
stress levels are below that normally considered critical 
for this cross-section. The maximum bearing stress was 0.55 
GPa (80 ksi) , maximum shear stress was 0.15 GPa (22 ksi) and 
the maximum tensile stress was 0.19 GPa (27 ksi) . 
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Close analysis of Figures IV-A-2-o(l-6) suggests bearing 
failure at the pin. Failure occurred at a bearing stress of 
65 ksl, which Is well below the value of 80 ksl that was used 
In design. This Is another Instance In which the ability to 
predict lug strength Is wanting, and Indicates that improved 
lug design Is an area which should receive additional study. 
This being the second 727 actuator attachment rib specimen 
to successfully sustain design load with adequate margin, 
we now consider this project concluded. 
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Figure IV-A-2-o(l). 

Failed Actuator: (a) Overall View, Side 1 
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Figure IV-A-2-o(2). 

Failed Actuator (b) , failure Site at Point E, View 1 



e IV-A-2-o(3) . 

Failure Site at Point E, View 2 
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Figure IV-A-2-o(4). 

Failed Actuator (d) , Failed Site at Point E, View 3 
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IV-B LOCKHEED L-lOIl ENGINE DRAG STRUT (CAPCOMP II) 

Senior Investigator: D. Goetachel 

1. Introduction 

Following a visit by Dr. R. Loewy to the Lockheed Cali** 
fornia Company on November 12, 1980 and subsequent meetings 
and discucsions in the NASA Headquarters office of the Pro- 
ject Officer, we jointly selected the drag strut of the wing- 
mounted, jot engine nacelles of the Lockheed L-1011 commer- 
cial transport aircraft for redesign in composites ao the 
second CAPCOMP project. A drawing of the strut, currently a 
steel forging is shown in Figure IV-B-1. The manner of its 
installation is shown in Figure IV-B-2, Both of these draw- 
ings, along with the first of the considerable information 
needed to perform this design, were graciously provided by 
the Lockheed Company. 


2. Status 

Preliminary studies were conducted to assess the magni- 
tude of the total task. Some of these considerations were 
reported in the progress report daced July, 1981. Other as- 
pects are suironarized here. 

Figure IV-B-3 is an extract of the design loads estab- 
lished by Lockheed for the part. Notice that compressive 
loads to 466.6 kips must be withstood. Accordingly, three 
basic compression failure criteria must be considered: 
material failure, Euler buckling and local buckling or 





Figure IV-B-2. 

Lockheed L-1011 Wing Pylon Structure Assembly 
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Figure IV-B-3. 

DESIGN FORCES FOR STRUT 



r 3 

DE SIGN ULTIMATE LOADS [iQ lbs.] 
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Condition 

Number 
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a 

P 

P 

z 

a 

P 

s 

m 



3 

-77.6 

-12.0 

-16.3 

79.3 

-118.68 

-19.69 

98.9^ 

21 

-456.6 

-0.7 

-96.2 

-466.6 

-6.92 

-236.76 

-229.84 

33 

341.2 

1.6 

71.9 

348.7 

15.82 

182.26 

166.44 

25 

-425.0 

3.8 

-89.5 

-434.32 

37.58 

-198.4 

-236.0 

34 

318.5 

-9.1 

67.1 

325.5 

90.0 

117.7 

207.7 

DESIGN LIMIT 

LOADS [ 10^ 







3 

-51.73 

-8.0 

-10.87 

52.87 

-79.12 

-13.13 

65.99 

21 

-304.40 

-0.47 

-64.13 

-311.07 

-4.61 

-157.84 

-153.23 

33 

227.4/ 

1.07 

47.93 

232.47 

10.55 

121.51 

110,96 

25 

-283.3 

2.5 

59.7 

-289.5 

25.1 

-132.3 

-157.3 

34 

212.3 

-6.1 

44.7 

217.0 

60.0 

78.5 

138.5 
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wrinkling. The formulas for handling the first two criteria 
are readily available in standard texts; however, formulas 
in a form suitable for design purposes for the local insta- 
bility of an anisotropic plate or shell are noc as readily 
available. One ot the first tasks for this project is to 
obtain or develop such equations for various boundary condi- 
tions. A search of the literature and design guides is now 
being conducted to see what formulas are available. If suit- 
able formulas are lacking, they will have to be developed for 
the pertinent and/or additional special cases of interest. 
Once these formulas are available, preliminary design concepts 
can be sized to determine their structural efficiency. Un- 
dergraduate student Frank DeTaranto, graduate student Shu- 
Ping Chen and postdoctoral fellow Christoph Muser have been 
pursuing this problem. So far both axisymmetric and nonsym- 
metric buckling formulas have been found for round tubes. 

3. Progress During Report Period 
Additional preliminary studies since the last reporting 
period have revealed the following facts. The size and load 
capacity of this part require special budgetary considera- 
tions. A full-sized strut will involve graphite/epoxy mate- 
rials costs of about $10,000; this severely limits the num- 
ber of design iterations that can be performed on full-sized 
parts. The Instron testing machine available at RPI has a 
maximum grip spacing of four feet and a load capacity of 
56,600 pounds. To test a full-sized part on-campus would, 
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therefore, require that a large and expensive mechanical 
lever mechanism be constructed. Testing off-campus would 
limit the number of Rensselaer personnel that could parti- 
cipate in the testing process. For these reasons we plan to 
perform most of the development on scale-model parts. Hope- 
fully then it would only be necessary to fabricate and test 
a single full-scale drag strut as a final demonstration of 
the design. A scale factor of one-quarter has been selected. 
This would make the scaled drag strut 2.5 feet long, which 
should fit in the Instron testing machine, even allowing for 
test fixtures. Material volume and costs will scale as the 
cube of the scale factor IC*)* * the material 

cost per strut to approximately $150. 

The remaining issue is whether the potential failure 
modes will all scale properly. Composite material failure 
stress is a composite material property; it should, there- 
fore, not change with scaling except for cure chemistry 
thermodynamics. For Euler buckling, only two geometric fac- 
tors - length and radius of gyration - are important, and 
they occur as a ratio of each other. Since they both scale 
with length, the Euler buckling stress should remain complete- 
ly unaffected by scaling. For local buckling of isotropic 
materials, the geometric factors also occur as ratios of 
length parameters, so this failure stress should a?.so remain 
unchanged by scaling for the anisotropic strut. With all the 
failure modes scaling in the same fashion, a full-sized strut 


should fail in the same fashion as the quarter-scale model, 
unless the ratio of fiber size to minimum part dimension 
should somehow become Important. 

If the failure stresses do not change with scaling, 
the failure load should scale with the cross-sectional area 
or the square of the scale factor. With a maximum axial load 
of 466.6 kips being scaled by (■x)’ , the maximum axial load 
for the scaled part should be 29.2 kips. T.hls is well within 
the capacity of RPI's Instron machine. The only remaining 
objection to sf:ale models would be difficulty in achieving 
minimum gauge thicknesses. It appears that the scaled mini- 
mum wall thicknesses will still allow for at least several 
plies of composite material, so this is not now expected to 
be a major problem. For the above reasons, quarter-scale 
model testing should provide realistic results, while allow- 
ing multiple interactions of the design-fabrication-testing 
cycle. 


4. Plans for Upcoming Period 
A critical aspect of the drag strut design is the way 
load transfer is effected at the connecting lugs at the ends 
of the strut. They are highly loaded, and in tension they 
have a very complex stress state. The lugs should be designed 
to fall within the same dimensional envelope as those on the 
production strut. Since these lugs are made from 200 ksi 
steel, meeting the geometric constraints could prove o be a 
difficult task even with a design making maximum use of 
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unidirectional graphite epoxy. We plan to attempt this first, 
but should it not be feasible, there are two other options. 
First, the lugs could be made from high strength steel and 
attached to a composite strut. Second, the attitude could 
be taken that this redesign is part of a redesign of wider 
scope that would allow composite lugs to be dimensioned as 
needed. 

The design of the lugs (see Figure IV-B-4) is important 
to the design of the rest of the strut, so this problem will 
be concentrated on first. Fixtures will be designed and con- 
structed for testing one side of a quarter-scale model lug. 
Wrapping th.“' graphite fibers around the mating pin would pro- 
vide the greatest tensile strength around the pin, but this 
approach has two drawbacks; there are problems with the outer 
fibers bunching up during curing and potentially very low 
transverse compressive strength next to the pins. Such prob- 
lems can be successfully dealt with, but they do compromise 
the effectiveness of this approach. Initially the suitabil- 
ity of a flat, symmetric angle-ply laminate will be investi- 
gated. Computer programs are being written to optimize the 
orientations and volume fractions of the plies. Scaled 
models of candidate lay-ups will then be tested to determine 
if this approach can be successful. 

It is noteworthy that three research areas of generic 
interest to composite structural elements have emerged from 
our Lockheed 10] 1 drag strut redesign efforts to date. They 
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are (a) the need for quantified information regarding St. 
Verant's principle for typical cross-sections of compression- 
tension-bending members using composite materials (semi- 
annual progress report, July, 1981); (b) the degree to which 
design development of major structural elements can be con- 
ducted with confidence at reduced scale; and (c) whether all 
composite, heavily loaded lugs, which will meet stringent 
over-all geometric envelope constraints, can be successfully 
designed and fabricated. We anticipate dealing with all of 
these important problems to some extent in the months ahead. 
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IV-C MECHANICAL JOINTS IN COMPOSITE S 

Senior InvestlgAtors: D. Goetschel 

R . Loewy 

The inunedlately preceding two sections of this report 
illustrate forceably the importance of being able to deal 
with mechanical joints in composite plates. In this section 
progress is reported on two continuing research projects in 
this area. The first involves optimizing the laminates 
around a hole - first open and later to be pin-loaded - 
using special variations of the material properties with 
distance from the hole. This is the work of post-doctoral 
fellow C. Muser. The second involves attempts to predict 
accurately stresses in the region of a pin-loaded hole in 
composite plates of typical, spatially uniform lay-ups. The 
second research is that of graduate student W. Kim. 

1. Analysis of Heavily Loaded Mechanical Joints 
a . Introduction 

Today's composite materials provide flexibility in 
choosing material properties for special design applications. 
In principle it is possible to have the fibers running in 
any direction, straight or curved, and with spatial varia- 
tions in their fiber directions. Structural elements made 
from these materials} can, therefore, be designed to specific 
properties at specific locations in order to increase load- 
carrying capacity-to-weiqht ratios. This aspect of fiber 
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reinforced plast-ics is not understood well enough to allow 
designing mechanical joints and openings in composite plates 
in a way which exploits the "tailored material" or'portuni- 
ties available with composites. A detailed investigation 
seems to be justified. 

b. Review of Previous Work 

A closed form solution was found by i>J. J. Hoff for the 
stress distribution in a circular plate with a central hole 
of radius r = 1 and an outer radius, r = R, under uniform 
uniaxial traction ,o, when the material is cylindrically or- 
thotropic and homogeneous. Figure IV-C-l-a shows the geome- 
try of such a plate. The maximum circumferential stress, 

at the boundary of the hole can be expressed as a func- 
tion of the outer radius, R, of the plate and the two mate- 
rial property parameters, and S. Here the first para- 
meter is defined as = ^rr'^^OO ' second parameter 

is S = -f . For any fixed 

gt?ometry of the plate (radius R fixed) , it is possible to 
draw lines of equal stress concentration-factors as a fvmc- 

ticn of the two paraneters S and S„ . (Stress concentration- 

rr 

factor is defined as the highest circumferential stress, 1 ^ 00 » 
at the edge of the hole divided by the average stress, 0 .) 

In other words, for a given geometry one can read the stress 
concentration-factor from a simple graph for any cylindrical- 
ly orthotropic material (See Figure 1-27, Semi-Annual Progress 
Report, July, 1981, pp. 66). 
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Part of the work in this area was to find the expressions 
for the stress distributions in the same kind of circular 
plate when the material properties are not constant but rather 
a function of the distance from the central hole. The con- 
stitutive equation used in this case is: 
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The two parameters and characterize the material 
properties at the edge of the hole and at the outer boundary 
of the disk. The parameter h is used to describe the change 
of material properties along the radius. The solution found 
is in the form of a power series of the radius ,r , and uses 
the exact closed-form solution for homogeneous materials as a 
starting point. The conditions for equilibrium and the 
boundary conditions are satisfied exactly. The condition for 
compatibility is satisfied U' the limit. It was found that 
the stated solution does not converge for every combination. 
Convergence limits were established. Several numerical re- 
sults were obtained. 


c. Progress During Report Period 

The first part of the current period was devoted to 
documenting and reporting the results obtained earlier. 

Progress was also made in applying the solutions found 
for the unloaded hole to the case of a circular plate pulled 
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at its outer edge r « R in one direction and pulled in the 
opposite direction by a bolt inserted in a central hole. 

The material properties are still o’lowed to vary along the 
radius^ r. The geometry usef* is shown in Figure IV~C-l-b. 

It is believed that cylindrical oithotropy - as in the 
case of the open hole - will result in low stress concentra- 
tions around the loaded hole as well. The aim of the work 
xs not to find the exact boundary conditions at the edge of 
the hole nor is it to find the proper failure criterion. 

The aim is rather to show the possible potential of cylindri- 
cal orthotropy as a means to design lugs more efficiently. 
Therefore, the calculations are kept strictly two-dimensional 
and the boundary conditions at the edge of the hole express 
only the presence of normal stress, o^, over one-half of the 
circumference. The common assumptions that the bolt fits 
snugly and that friction between the bolt and the plate does 
not exist were also made here. The boundary conditions can 
be expanded iiico a trigonometric series. It is possible to 
adjust the stress function found for the empty-hole case so 
that it can be used to solve the present problem. Moreover, 
it seems possible to solve for any boundary conditions which 
can be expressed in trigonometric series, so long as con- 
vergt-nce is provided. As of this writing, no numerical re- 
sults have been obtained. 
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d. Plans for Upcoming Period 

In the months ahead we plan to continue the investiga- 
tion of stress distributions in circular, bolt-loaded plates, 
the material properties of which are cylindrically orthotro- 
pic and vary along the radius. In addition, we expect to 
make presentations of and to publish the results found for 
the open-hole case with radially varying material properties. 

e. Current Publications or Presentations on this Subject 
by Professor Hoffs 

"Stress Concentrations in Cylindrically Orthotropic 

Composite Plates with a Circular Hole" 

Published in the Journal of Applied Mechanics, 
September, 1981, Vol. 48, pp. 563-569. 

Presented at the Joint ASME/ASCE Mechanics 
Conference, Boulder, Colorado, June 23, 1981. 

by Dr. Muser: 

"Cylindrically Orthotropic Circular Plates with a 

Concentric Hole and with Properties Varying Radially" 

Thesis for the degree. Doctor of Engineering, 
Rensselaer Polytechnic Institute, August, 1981. 

2. Pin-Loaded Holes in Uniform Composite Plates 
a. Introduction 

To optimize a mechanical joint in composites, it is nec- 
cessary to be able to predict the stress dirtributions around 
the hole. In this approach to solving this problem, analysis 
has been restricted to a single hole, plane problem, as a 
necessary first step. 

Even this simplified problem falls in the category of 
mixed boundary value problems with unknown boundary conditions 
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on a part of a bounda :y. The amount of contact between the 
loading pin (assumed rigid) and the hole as it deforms under 
load is thus unknown. By introducing a semi-empirical analy- 
sis method, which we called the "mixed method", as explained 
in previous reports, the boundary conditions can be deter- 
mined. 

One of the most important results of this reporting per- 
iod is that the possibility has emerged that resolution of 
the contact problem may be obtained - without the "mixed 
method" - by using interactive graphic computer terminals. 

This approach, as will be explained in the subsequent section, 
is based on the idea that an originally smooth curve, in or 
on the boundary of the material, must remain smooth after de- 
formation, provided that the surface applying the load is 
also smooth. 

b. Status 

As noted in earlier reports, once the extent of the re- 
gion of contact between the loading pin and the loaded hole 
is known (originally undetermined for clearance fit cases) , 
the stresses in the domain can be obtained by various means. 
For exeunple, a finite element analysis can be performed with 
displacement boundary conditions for points located in the 
contact region. These boundary conditions are determined by 
the rigid surface applying the load. During the previous 
period, a system of linear constraint equations was applied 
to the nodal points in the presumed contact region for this 
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purpose. The correct extent of the contact region (defined 
by the contact angle 0 ) was inferred by comparing the maxi- 
mum strain contours from several such finite element analy- 
ses with different values of 0 with the isochromatic fringes 
from corresponding photoelastic experiments. Where the maxi- 
mum shear strain distributions were well-matched, the assumed 

value for 0 , as well as the corresponding stress solutions, 
c 

were taken as also correct. 

The feasibility of this method for this particular prob- 
lem has been discussed in previous reports. We note here, 
however, that this "mixed method" can be expanded for more 
complicated problem situations; frr example, multiple pin 
connections; and second, that tb® results are not subject to 
uncertainties as to the adequacy of the finite element model 
or other assumptions, since the method is self-checking pre- 
cisely because of its semi-empiricism. 

c. Progress During Report Period 

The linear constraint equation used to prescribe the 
motion of a node point on the boundary within 0^ has been 
changed during this report period. It was derived as fol- 
lows . 

Let the undeflected position of a node point, p, be x 
and y in the cartesian coordinate system centered on the hole 
(Figure IV-C-2-a) . If p is within 0^, for a particular level 
of applied load, then the final location of p after deforma- 
tion is somewhere on the surface of the rigid, fixed pin; 



125 


i.e., at a point described by 

(x + u) * + ((y - 6) + vl* » r* 


( 1 ) 


where, u and v are displacements of p, in the x and y direc 
tions, respectively, 6 is the clearance and r is the radius 
of the pin. Equation (1) can be rewritten by expanding and 
neglecting u^ and v^, since they are small, first order. 

X* + 2ux + (y - 6) * + 2v(y - 6) “ r^ (1') 


Dividing both sides of Equation (1^) by Jtk} + (y - 6) ^ , 
noticing that 


/x* + (y - 6) 


* cos 0 


and 


sin e 


/x’ + (y - 8) 


and rearranging, one obtains 


u cosO + V sin 0 = 


- [x^ (y - 6) M 

2/x* 4 (y - 6) ^ 


( 2 ) 


This equation is the form in which the constraint equettion 
was used. 

If 6 = 0 in Equation (2) (i.e., no clearance), then 
Figure IV-C-2-a makes it clear that the right-hand side of 
(2) is zero, so that 


u cos 0 + V sin 0 « 0 


(3) 
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This is the constraint equation for the case of no clearance 
and is consistent with constraints used in previous reports. 

Also during the reporting period, two photoelastic spe- 
cimens were manufactured and tested. The overall dimensions 
of the new specimens are almost identical to those of previ- 
ous specimens, except for the thickness. As shown in Table 
IV-C-2-a, the total number of layers was Increased from 18, 
as used in the old specimens, to 23 for B-type lay-ups (i.e., 
(±45“/0®l ). For this particular laminate, n was intended 

to be 8, since the design called Tor 24 layers, but one 0® 
layer was inadvertently left out during the lay-up process. 
However, symmetry was maintained, since it was a center 0® 
layer that was missing. The resulting thicknesses of plates 
after curing were 0.132 inch for l±30®/0®]gg and 0.127 inch 
for B-type plates. This matched well with what would be ex- 
pected from manufacturer's standards for graphite/epoxy pre- 
preg tapes, 

Plate thickness increases were made the subject of in- 
vestigation since earlier results suggested that outer layer 
delamination and subsequent buckling of separate fibers were 
mechanisms in bearing failures. This effect may be made a 
smaller fraction of total strength as the total number of 
layers increases. 

One may see from Table IV-C-2-a that the bearing strength 
of a composite plate is strongly dependent on the total thick- 
ness but not so much on the ratio of the hole diameter to the 
thickness (D/t) . More precisely, the D/t of Experiment II 
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TABLE IV-C-2-a 

STRENGTH AS INFLUENCED BY D/t FACTOR, (±45®/0‘’l 


Experiment 



Diameter 

D 

(inch) 


0.25 


.5 


0.4995 


Number 

of 

Layers 


18 



0.103 


0.103 




D/t 

Average 
Bearing Stress 

(Psi) 

2.43 

71,706 

4.85 

70,526 

3.93 

94,488 



0.4995 


0.127 


0.9965 0.5 


0.132 
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(4.85) is twice as big as that of Experiment I (2.43) > but 
the average bearing stresses at failure are roughly the same 
for both specimens and their thicicnesses are exactly the 
same. On the other hand, the roughly 30% thickness increase 
of Experiment III resulted in more than a 50% increase of 
failure load compared to Experiments I and II, although the 
Experiment III plate's D/t fell between those of the earlier 
plates . 

While this result indicates that the bearing strength 
of a composite joint depends highly on the overall thickness 
of a plate instead of on D/t ratio, we believe that there 
will be a limit to the total number of layers alaove which 
increasing the number of layers will affect bearing strength 
only negligibly. This follows since the loss of strength of 
the more easily delaminated outer layers becomes a smaller 
percentage of the total strength of the full number of lay- 
ers as the laminate thickness grows, as noted earlier. 

This has important meaning from a model scaling point 
of view. Apparently, one must take this effect into account 
whenever scaling up (or down) pin-loaded joints in composites. 

Other progress during the subject period includes im- 
provement in the manufacturing processes of photoelastic spe- 
cimens. Grid nets were engraved on the back side of the 
photo-elastic coatings, prior to bonding on the test speci- 
mens, for easy fringe position determination. Special atten- 
tion was given to specimen dimensions to insure symmetry. 


Final measurements showed that the eccentricity of the loca- 
tion of the hole was less than 0.001 inch for both specimens. 
Hole diameter was checked in several directions usiiKi a flex- 
ible radius gauge > which^ when used skillfully to check round- 
nes8» can detect 0.0005 inch differences. The differences 
between maximum and minimum diameters in the holes in our 
specimens were too small to be detected using this instrument. 

The loading pin used in these tests was manufactured so 
that the diametral clearance would be 0.0015 inch after all 
the necessary corrections and measurements had been made. 

Due to the increase of laminate thickness, several metal 
parts of the teat fixture had to be heat treated to give 
higher strength (see Figure IV-C-2-b) . 

Photoelastic experiments were performed on an Instron 
testing machine. At evr«ry additional 500 pounds of load, 
isochromatic pictures were taken until the specimen failed. 

The failure loads were 6,000 pounds for B-type specimens and 
5,860 pounds for F-type (lt30yo"lgg) laminates. It is inter- 
esting that the total number of layers of Specimen B was one 
layer less than that of F-type, but still showed higher 
strength. Table IV-C-2-b shows the detailed dimensions meas- 
ured for both specimens. 

Additional finite element analyses have been performed 
in this period, following the same process that has been de- 
scribed in several previous reports. Figure IV-C-2-c shows 
the isochromatic fringes of Specimen B at 4,000 pounds, and 
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Figure IV-C-2-d shows the maximuin shear strains resulting 
from the finite element analysis. The value of the maximum 
shear strains for corresponding contours are given at the 
upper right-hand side of Figure IV-C-2-d. These values a"' 
respond to those of Figure IV-C-2-c and will oe seen to check 
very well. Figure IV-C-2-e is a scaled-up picture of Figure 
IV-C-2-d. The maximum shear strain distributed in both are 
almodt identical to those in Figure IV-C-2-f, except for 
minor details. This leads to the conclusion that the con- 
tact angle of tnis specimen with 0.0015 inch no-load clear- 
ance between the loading pin and hole is about 75 degrees at 
a load level of 8,000 pounds. 

The same procedures have also been completed for load 
levels of 1,000, 2,000, 3,000 and 5,000 pounds. The results 
allow plotting the relationship between applied load and 
contact angle as shown in th« Figure IV-C-2-f. By extrapo- 
lating and interpolating this curve, one can obtain the con- 
tact angxe at any load level including the failure point. 

This enables us to examine the stress situations at the maxi- 
mum load with considerable confidence. 

The ranges of contact angles shown in Figure IV-C-2-f 
are those within which the contour plots of maximum shear 
strain are too similar to allow selecting one as being a bet- 
ter match thiin the others. However, once those ranges of 
contact angle are put on one plot, as in Figure IV-C-2-f, 
and assuming the correct value is closer to the middle of 
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of the range than the extremes, one can easily choose the 
best contact angle f ' r analysis. 

In the process of performing these experiments and cal- 
culations, it has become clear that, even beyond satisfying 
boundary constraint conditions on displacements, curves ori- 
ginally smooth in the material domain must remain smooth 
after elastic deformation. Thus, it may be possible to find 
a correc' contact region by plotting deformed shapes and ex- 
amining the smoothness of their contours. 

One aspect is encouraging even at this point; the lowest 
maximum shear strain values were observed when the deformed 
shape was best matched to such a smoothness condition. 

d. Plans for Upcoming Period 

The "mixed method" has proven to this point to be a 
practical and reasonably good research approach to a compli- 
cated and uncertain boundary value problem. The alternative 
smoothness condition mentioned must be examined further. If 
it shows itself to be reliable, then it will be possible to 
produce stress data for various combinations of fiber arrange- 
ments without the time and labor of photoelastic experiments. 
We will investigate this and, if successful, try to apply a 
failure criterion for pin-loaded composite jcints in beating. 
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IV-D COMPUTER SOFTWARE DEVELOPMENTS 

Senior Investigator: M. S. Shephard 

1. Introduction and Status 

The objective of this portion of the project is to pro- 
vide advanced and specialized computer "tools" for the anal- 
ysis and design of composite structural elements. The major 
thrust to date has been and continues to be in the finite 
element area, with effort directed at implementing finite 
element analysis capabilities and developing interactive 
graphics preprocessing and postprocessing capabilities. Re- 
cent efforts have been directed at completing the two-dimen- 

* 

sional interactive graphics programs in the POPES system 
and initial development of a three-dimensional program. 

2. Progress During Report Period 

Effort during the current reporting period on POPES 
software has concentrated on the following: 

(a) Development of an edge load editor for the 
two-dimensional attribute editor. 

(b) Initial development of a surface load edi- 
tor for the two-dimensional attribute edi- 
tor . 

(c) Initial development of a preprocessor for 

** 

general three-dimensional shell surfaces 


POPES (People Oriented Finite Element Software) is RPI's 
finite element software system. 

Funding for this work is external to the composites pro- 
grctm and through the Graphics Center's Industrial Asso- 
ciates Group. 
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The function of the edge load editor la to allow the 
user an efficient means to define edge tractions. This seg~ 
ment of the attribute editor has features that allow the user 
means to define his edge loads In terms of components normal 
to an edge, tangent to an edge or at any given angle to an 
edge. The load can be uniform along an edge or have a lin- 
ear, quadratic or cubic variation. All such loads, however, 
are understood to lie In the plane of the object. Figure 
IV-D-1 has three edge loads applied. The upper left edge 
load Is a quadratlcally varying normal edge load while the 
second Is a constant load acting at an angle to the edge. 

The length of the arrows used to display the load Indicates 
the relative magnitude of the load, and the load direction 
is indicated by the orientation of the arrows. The third 
load is a linearly varying tangential load on the horizontal 
line at the top of region 2. The relative magnitude of tan- 
gental loads is indicated by the size of the arrow heads 

* 

used to display it. Figure IV-D-2 shows the use of the ECHO 
function to check the edge load. 

The function of the surface load editor is to allow the 
user an efficient means to define loads acting over and nor- 
mal to the surface of the object. Figure IV-D-3 shows the 
surface load menu page with a constant normal surface load 
applied to the right-hand region. In Figure IV-D-4 a linearly 
varying load has been applied to the next region, which 
matches the constant load on the neighboring region. Figure 
1 

ECHO is a function by which the operator can call back for 
review any part of data or program. 
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Xv-D-5 ahow8 a complete set of surface loads on the same 
structure. The user can also check the values of the surface 
loads applied using the ECHO function as shown in Figure IV** 
D-6. 

A surface mesh preprocessor for the generation of fi- 
nite element meshes for three-dimensional shell structures 
is being developed. As with the other graphics progreuns in 
the POFES system, the objective of the surface mesh prepro- 
cessors is to offer the user an interactive graphic tool 
that is easy to use. 

The process of defining a finite element surface mesh 
begins with defining the space the user wishes to work in 
and defining the three-dimensional boundary curves. Several 
features have been built into the program to aid the user in 
the input of the three-dimensional boundary curves. The most 
important of these gives the user a means to obtain the exact 
values of specific three-dimensional locations by "Hit-test- 
ing"* against the two-dimensional projection of various 
points displayed on the screen. The points against which the 
•'hit-test'' may be conducted include the ends of previously 
defined curves, the key points used to define the curves and 
the currently active grid plane which is user-selected from 
the user-defined three-dimensional grid. The user is also 
given three different methods for rotating and translating the 


"Hit-testing" is a process in which the coordinate values 
given by the interactive input device are set to those of 
an existing location, if the "distance" between the two 
locations is less than some prescribed value. 
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object, along with the ability to zoom in or out, to trans- 
late and change the scale of the display. Figure IV-D-7 
shows a set of cur^/es interactively input using the feature 
just defined. The user can c’toose to use either cubic 
splines or straight lines to define the boundary curves. 

The curves are then set to the mesh-generating portion 
of the progreun. Here the user Interactively breaks up the 
curves into the desired number of element edges, specifies 
the desired element type and selects the curves that make 
up the surface mesh. Figure IV-D-8 shows the surface mesh 
generated using the boundary curve shown in Figure IV-D-7 . 
Figures IV-D-9 and 10 show two views cf another surface mesh 
example. 


3. Plans for Upcoming Period 
During the next reporting period the final programming 
on the attribute editor will be finished and complete user 
documentation prepared. Other work on POFES software will 
include further development of the surface meshing program, 
along with improvements to the inhouse general purpose fi- 
nite element analysis program to be carried out by a small 
group of students working on course projects. 

4 . Current Publications or Pre sentations by 
Professor Shephard on this~Subject 

1) "The University's Role in Meeting Industry's CADCAM Needs" 

Published in "Education and Industry - A Joint En- 
deavor", Vol. 1, Proc. ASEE Nat. Meeting, June, 

1981, pp. 232-235. 
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Figure IV-D—7. Definition of Three-Dimensional Boundary Curves 
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Figure IV-D-9. Surface Mesh Example 
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2) ”A Osncraliicd Two-Dlmanslonal, Graphical Plnita Blamant 
Praprocaasor Utilizing Diacrata Transfinita Mappings*', 
(with R. B. Habar, J. P. Ahal, R. M. Gallaghar and D. P. 
Graanbarg) . 

Publiahad in Int. J. Num. Math. Bngng., Vol. 17, 

July 1981, pp. 1015-1044. 

3) "Pinita Blamant Modal Ganaration" 

Praaantad at Battalia Columbua Laboratories, Ohio, 

Saptambar, 1981. 

4) "Pinita Element Solutions to Point-Load Problems", (with 
R. H. Gallaghar and J. P. Abel). 

To be published in J. Eng. Mach. Div., ASCE, Vol. 

107, No. EM5, October, 1981, pp. 839-850. 

5) "The Finite Element Modeling Process - Will It be Auto- 
mated" 

To be published in "New and Puture Developments in 

Commercial Finite Element Methods", J. Robinson, 

Editor, October, 1981, pp. 451-468. 

6) "Adaptive Finite Element Analysis Accounting for Material 
Interfaces", (with S. M. Yunus). 

To be presented at the ASCE Conference, St. Louis, 

Missouri, October, 1981. 
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V-A INITIAL SAILPLANE PROJECT: THE RP-1 


Senior Investigator 8 : 


P. Bundy 

R. J. Diefendorf 

H. Hagerup 


Flight tests of the first ultra-light all-composite 
glider resumed in the Summer of 1981, with the number of 
flights of the RP-1 reaching a tox,t.l of thirty-six at thf- 
end of the season, h camera was installed to record flight 
instrument readings, and several flights were conducted in 
the still air of the early morning hours, so that a glide- 
polar could be determined. The preliminary results are 
shown in Figure V-A-1, from which the 1981 configuration of 
the RP-1 can bo seen to have a best glide ratio of 11.3 : 1 
and a minimum rate of sink of 4.2 fps. 

It is necessary to bear in mind that during these 
flights the aircraft was equipped both with an unfaired 
landing wheel and with a camera mounted in an exposed posi- 
tion. It is estimated that the cleaner, 1980 configuration 
had a parasite drag coefficient about fifteen per cent lower. 
The aircraft has also become seventy pounds heavier since 
then. Consequently, an effort will be made during this win- 
ter to reverse some of these adverse aerodynamic developments 
by installing the necessary fairings under and behind the 
pilot. Not the least of these aerodynamic improvements will 
be a closure for the open fuselage bottom. This opening was 
originally to allow foot-launching, but the degree of success 
obtained to date with winch-launching has made that mode of 



airspeed ( mph ) 



Glide Polar - Preliainary Data 
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launching unnecessary. Accordingly, the fuselage closure Is 
possible, and this should result in much lower induced drag 
and less lift- loss associated with "leakage” inboard of the 
wing roots. We expect to be able to report a best glide 
ratio close to fifteen for the RP-1 by the end of the Summer 
of 1982 and fairings which add very little to the aircraft's 
gross weight. 
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V-B SECOND SAILPLANE PROJECT i THE RP-2 

S«nlor Investigators: P. Bundy 

R. J. Diefendorf 
H. Hagerup 

This project encompasses the design, fabrication and 
testing of an ultra-light, high performance, all composite 
sailplane. A number of sophisticated design aspects have 
been selected to challenge Innovation in fabrication pro- 
cesses for low cost aerospace structures beyond tho level 
encountered in the first project, which led to the RP-1. The 
RP-2 design is now complete and includes such features as a 
partial span split-flap, horizontal tail with field splices 
at roughly 2/3 span, double-taper wing planform and tail 
length to mean aerodynamic chord ratio of 5.3. 

Fabrication has progressed to the point where empennage 
and wings are now ready for static proof testing to design 
loads. Since the last report, the partially completed frame- 
work of the wings was inspected by the FAA, and the quality 
of the work was approved. 

The final design for the L-section stringers, which 
stiffen wing skins against buckling, incorporates a hybrid 
lay-up with the following stacking sequence: 0** graphite, 

±45** Kevlar cloth and 0^ graphite. These stringers are suf- 
ficiently stiff to eliminate the use of additional shear 
webs between the spar and the leading edge, looking forward, 
or between the spar and the flap/aileron hinge line, looking 
aft. The low linear density of 23 g/m achieved for these 
stringers results in a significant weight saving. 
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The wing skins were Installed after FAA inspection, 
without the use of external molds of any kind during wing 
construction. Tolerances of the finished wings were deter- 
mined to be completely satisfactory, confirming our expecta- 
tions of the labor-saving wing construction technique used 
for the first time on the RP-2. The all-graphite split 
flaps have been designed with skins lald-up as 4-plies, 

(±45**) , and assembled on triangular ribs with foeun spacers 

s 

to maintain cross-sectional shape. Bell cranks, pushrodi 
and torsion-rods for aileron and flap controls have been in- 
stalled, with removable, covered inspection ports for ease 
of access at the appropriate locations. 

The most interesting structural feature of the fuselage 
design for the RP-2 is the use of two large, curved hybrid 
graph:. te/Kevlar foam sandwich beams which carry the loads 
from towline, pilot weight and landing wheel reactions into 
the center box which connects wings and tailboom. One such 
full-scale prototype curved beeun, referred to as a "tusk" 
and designed to carry 500 pounds with a safety factor of 
three, was tested to destruction with very satisfactory re- 
sults; failing in shear at 1,650 pounds. The production 
becuns, consequently, were laid-up precisely as the prototype 
and have now been installed, thus completing the structural 
framework of the aircraft. Control linkages within the 
fuselage have also been essentially completed, and the only 
major fabricating task remaining on the RP-2 is the construction 


of the fuselage fairing enclosing the pilot and framework. 
This task is expected to be completed during the Fall and 
Winter of the next reporting period, in anticipation of 
initial flight of the RP-2 in the Summer of 1982. 
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PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy, Robert G.» Ph.D. 
Wiberley, Stephen E.» Ph.D. 

Senior Investigators 

Brunelle, E. J. , Jr., SC.D. 
(Aeroelastic and structural 
design and analysis, applied 
mechanics of composite 
structures, research)* 

Bundy, F. , Ph.D. 

(Physical chemistry and 
CAPGLIDE) * 

Diefendorf^, R. J., Ph.D. 
(Fabrication, CAPGLIDE, 
fiber behavior, research)* 

Feeser^, L. J., Ph.D. 
(Computer applications and 
graphics, computer aided 
design, optimization)* 

Goetschel, D. B., Ph.D. 
(Structural analysis and 
design, CAPCOMP ) * 

Hagerup, li. J., Ph.D. 
(Aerodynamics, configura- 
tion, pilot accommodation, 
CAPGLIDE) * 

Hoff, N. J., Ph.D. 
(Structural design and 
analysis, CAPGLIDE)* 

Krempl, E., Dr.Ing. 

(Fatigue studies, 
research) * 

Scarton, H., Ph.D. 

(Acoustic emission NDE) * 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 


Associate Professor of 
Aeronautical Engineering 


Research Professor of Materials 
Engineering 

Professor of Materials 
Engineering 

Professor of Civil Engineering 


Assistant Professor of 
Mechanical Engineering 

Associate Professor of 
Aeronautical Engineering 


John A. Clar)c and Edward T 
Crossan Professor of Engineer- 
ing 

Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Associate Professor of Mechan- 
ical Engineering and Mechanics 


Fields of Speciality. 

^Member of Budget Committee together with co-principal investi- 
gators. 



158 


Senior Inveatiqatora 

Shephard, M. , Ph.D. 
(Computer graphics, finite 
element methods) * 

Sternstein^, S. S., Ph.D. 
(Failure analysis, matrix 
behavior, research)* 


Assistant Professor of Civil 
Engineering 

William Weightman Wallcer 
Professor of Polymer Engineer 
ing 


Research Staff 

Manager & Master Technician, Composites Laboratory 
Paedelt, Volker 


Graduate Assistants 


Bobal, Gail, M.E. 

Chang, Chi-Min, M.S. 

Chen, Chikuang, B.S. 

Chen, Shu-ping, B.S. 
Cummings, Delwyn, B.S. 

Das, Nripendra, B.Tech. 
Helmer, James, B.S. 

Kim, Wonsub, M.S. 

Lumban Tobing, Frida, M.S. 


Muser, Christoph, Dipl .Mech. Ing 
Niu, Tyan-Min, M.S. 

Oyibo, Gabriel, M.S. 

Uzoh, Cyprian, B.S. 

Valicenti, Raymond, B.S. 
Winckler, Steven, M.E. 

Yang, Phillip, M.S. 


Undergraduate Assistants - Seniors 


Adel, Paul 
Applewhite, Keith 
Bertolazzi, Andrew 
Blake, William 
Bristol, Brian 
Cox , Mary 
DeTaranto, Francis 
Dhimitri, Willieun 


Fairchild, Kendall 
Fisher, Steven 
Niederer, Melvin 
Northrop, Steven 
Putnam, David 
Staniorski, Anthony 
Vaney, Philippe 


Fields of Speciality. 

^Member of Budget Committee together with co-principal invest! 
gators. 




Undergraduate Aesietants - Juniors 


Baucoro^ Allan 
Carbone « James 
DeMint, Thomas 
Ehrgott, Darhl 
Ficarra, Robert 
Kirshenbaum^ Susan 


Kuntsmann, Debra 
Lagrenlexe, Suzanne 
Lundquist, Cheryl 
Malusa, Stephen 
Melendez, Hermenegll^o 


Undergraduate Assistants - Sophomore 
Schroeder, Hans 


PRECEOmO P»OE BLANK NOT «« 


161 


AUTHOR INDEX 


Brunelle, E. J., Jr. 

Bundy , F 

Diefendorf, R. J. .. 
0OG^8Ct^6lf D. B. . . . 

Uagerup, H. J 

Krempl, E 

DoGwy f R. G. 

Scar ton, H. A 

Shephard, M. S 

Sternstein, S. S. .. 


Page (a) 
111 

149 , 152 


21 , 

149 , 

152 

73 , 

105 , 

115 


149 , 152 
47 

73 , 115 
73 

58 , 64 , 137 
49 







